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NOMENCLATURE 

Material  Properties 
D  Grain  width 

Eg  Elastic  modulus 

E{j  Contribution  to  plastic  modulus  in  the  grain  from  slip 

Epg  Plastic  modulus  in  the  grain 

Ep^  Plastic  modulus  in  the  matrix  surface  layer 

Eq&Eq'  Material  parameters  used  to  express  in  terms  of  matrix  and  slip 

components 

Y  Crack  growth  rate  exponent 

Yp  Y2  Material  parameters  used  to  describe  slip  band  activation 

Tr3»  Y4  Material  parameters  used  to  describe  slip  band  hardening 

Oy  Bulk  0.2%  yield  strength 

Lower  flow  stress 
Upper  flow  stress 


Constitutive  Parameters 
Og  Constraint  in  the  surface  layer 

Qg  Constraint  in  the  grain 

Local  plastic  strain  in  the  grain 
Local  plastic  strain  range  in  the  grain 
Acp  Local  plastic  strain  range  in  the  matrix  surface  layer 
Acp  Plastic  strain  range  for  local  initiation  of  slip  bands 
Cj.  Measured  residual  strain  over  one  cycle 

Og  Local  stress  in  the  grain 
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Aog  Change  in  Icxral  stress  in  the  grain 
o®  Applied  stress 

°*max  Maximum  applied  stress 
o^min  Minimum  applied  stress 
p  Linear  density  of  active  slip  bands 

N  Fatigue  cycles 

n  Number  of  cycles  to  produce  slip  bands 

n^  Number  of  cycles  in  stress  spectra  (no  plastic  deformation) 

np  Number  of  cycles  in  plastic  strain  spectra 

n  Number  of  band  nucleation  sites/area  of  each  type 

Damage  Parameters 

a  Crack  depth 

Aa  Incremental  crack  depth  growth 

A6  Crack  tip  opening  displacement 

aH  Incremental  change  in  band  hardening  level 

Al  Incremental  change  in  band  nucleation  damage 

m'  Number  of  blocks  to  nucleate  a  slip  band 

P(  )  Probability  distribution  of  damage  parameters 

c  Linear  accumulation  of  damage 

aK  Stress  intensity  range 

Nomenclature  in  Papers 

a  Constraint  factor 

Ep  Plastic  modulus 

Eg  Elastic  modulus 
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Acp  Increment  in  plastic  strain 

a  Applied  stress 

Ao  Change  in  local  stress 

°mean  ^^®an  stress 

Wq  Width  of  applied  stress-local  strain  hysteresis  loop  at  zero  load 


3 

C8577D/bje 


NADC-89044-60 


Rockwell  International 

Science  Center 


SC5470.fr 


1.0  INTRODUCTION 

Conventional  cumulative  damage  calculations  based  on  macroscopic  strain 

range  work  reasonably  well  for  predicting  the  effect  of  load  spectra  on  the  fatigue 

lifetime  of  some  alloys,  especially  steels. These  same  procedures  are  notoriously 

inadequate  for  many  other  structural  materials,  including  aluminum  alloys.  For 

aluminums,  the  reason  the  conventional  analysis  fails  is  known.  The  deformation  that 

leads  to  cracking  occurs  in  localized  areas,  i.e.,  individual  grains,  at  the  alloy 
U  5 

surface.  ’  Hence,  the  local  stresses  and  strains  differ  substantially  from  the  applied  by 
amounts  which  depend  on  load  sequence. 

The  purpose  of  the  research  under  this  one-year  contract  was  to  advance  the 
understanding  of  the  deformation  processes  that  control  the  local  stresses  and  strains  in 
order  to  formulate  a  rational  model  of  crack  initiation  under  spectrum  loads.  Experi¬ 
mental  techniques  of  high  spatial  resolution  strain  measurement  developed  at  the  Science 
Center  were  used  to  investigate  the  mechanica'  behavior  of  individual  constrained  grains 
in  A1  2219-T851.  Local  mechanical  properties  and  yield  criteria  can  be  deduced  from 
measurements  of  the  local  stress-strain  hysteresis  loop  width.  Some  of  our  recent 
findings  have  changed  the  structure  of  the  model  initially  envisaged  for  the  A1  2219-T851 
material  studied: 

1.  Fatigue-induced  lowering  of  the  local  yield  strength  in  isolated  micro¬ 
plastic  grains  in  A1  2219-T851  is  not  the  result  of  a  change  in  flow 
stress.  Instead,  the  plastic  modulus  of  the  grains  decreases  progressively 
with  fatigue.^  A  distribution  in  the  strain  range  required  for  slip  band 
nucleation  is  needed  to  explain  the  load  history/stress  amplitude 
dependence  of  grain  softening. 

2.  The  decrease  in  stress/strain  loop  width  seen  in  grains  later  in  lifetime 

results  from  strain  hardening.  Further,  if  individual  grains  are 

completely  hardened  at  one  cyclic  stress,  they  will  revert  to  a  softened 
state  if  the  cyclic  stress  is  raised.  The  hardening  rate  is  faster  in  grains 
in  which  the  local  plastic  strain  is  large.  Hardening  is  probably  not  due 
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to  cross  slip,  but  rather  to  cumulative  damage  in  each  slip  band  that 
ultimately  deactivates  slip  in  the  band. 

3.  The  two  stages  of  flow  found  in  A1  2219-1851  involve  different  constitu- 

4  5  7 

tive  relationships.  ’  ’  The  high  strain  hardening  flow  at  small  stresses 

b 

is  kinematic  and  we  believe  this  occurs  in  A1  2219-T851  in  a  thin  layer 
at  the  surface.  A  kinematic  behavior  could  imply  that  there  exists  a 
short  range  internal  constraint  of  the  deformation  which  creates  a  local 
back  stress. 


4.  Flow  above  the  upper  flow  stress  in  A1  2219-T851  is  stationary,  sugges¬ 
tive  of  unimpeded  deformation  until  the  slip  bands  reach  the  grain 
boundaries.  A  model  of  the  stress-strain  behavior  of  the  two  flow  stress 
systems  has  been  proposed  by  viewing  the  material  as  a  microscopic 
composite.^ 

5.  Consistent  with  the  picture  that  softening  involves  the  creation  of  slip 
bands  in  a  surface  layer,  which  then  penetrate  the  grain  body,  is  the 
finding  that  all  A1  2219-T851  grains  soften,  but  that  the  initial  strains  in 
large  grains  are  much  larger  than  in  smaller  grains. 

6.  A  convenient  engineering  definition  of  crack  initiation  in  A1  22i9-T851  is 
shown  to  be  the  cycles  required  for  the  localized  plastic  flow  in  a  grain 
to  assist  microcrack  propagation  to  the  first  subsurface  grain  boundary. 


In  the  next  section,  we  discuss  these  findings.  We  then  proceed  to  a  descrip¬ 
tion  of  a  lifetime  initiation  model  which  supposes  a  microscopic  inhomogeneity  of  the 
deformation  of  individual  A1  2219-T851  grains  to  obtain  expressions  to  calculate  the 
locally  evolving  mechanical  properties  and  the  resulting  crack  initiation  rates.  This 
model  is  phenomenological.  It  replaces  several  areas  of  analysis  where  mechanistic 
insight  is  still  lacking  with  empirical  expressions  that  can  now  be  directly  calibrated  with 
microscopically  measured  deformation  parameters.  Perhaps  in  the  future,  calibration 
can  be  indirect  through  lifetime  response  to  standardized  load  sequences. 
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2.0  OBSERVATIONS  ON  LOCAL  DEFORMATION  IN  AI  2219-T851 

Questions  raised  in  our  initial  proposal  were  the  following: 

1.  Is  the  slip  band  model  we  had  proposed  (under  NSF  contract)  for  grain 
softening  correct? 

Results  in  the  paper  presented  in  Section  2.1  indicate  that  some  of  the 
essential  features  of  the  deformation  are  correctly  predicted.  For 
example,  softening  is  associated  with  a  decreasing  plastic  modulus,  as 
expected  when  the  slip  band  density  increases.  The  fine  details  are 
different  from  that  originally  imagined,  however.  For  instance, 
deformation  below  the  upper  flow  stress  strain  hardens  with  fatigue, 
possibly  because  the  slip  bands  reduce  the  volume  of  the  surface  layer 
which  can  deform.  To  describe  the  softening  rate  vs  local  strain  range, 
we  propose  an  empirical  model  (Section  2.3)  requiring  calibration  by 
microscopic  deformation  measurements. 

2.  Does  deformation  associated  with  each  of  the  two  separate  flow  stresses 
in  Al  2219-T851  obey  the  same  yield  locus  behavior? 

No.  As  described  in  the  paper  presented  in  Section  2.2,  deformation 
between  the  lower  flow  stress,  and  o^,  the  upper  flow  stress,  is 
kinematic.  Above  the  deformation  flow  criterion  is  stationary.  This 
is  a  further  confirmation  that  flow  below  is  localized  in  the  grain 
matrix,  and  above  occurs  in  long  ranged  bands  which  reach  the  grain 
boundaries. 

3.  Can  we  ignore  the  local  hardening? 

It  could  be  seriously  over-conservative  to  do  so.  For  high  stress  ampli¬ 
tudes  near  yield,  correction  for  hardening  cannot  matter  because  crack 
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initiation  occurs  in  the  large  grains  before  they  have  time  to  harden. 
But,  at  lower  amplitudes,  the  much  slower  rate  of  hardening  in  some 
smaller  grains  makes  them  the  potential  cracking  targets,  in  comparison 
to  the  very  large  grains  which  are  protected  from  cracking  by 
hardening.  (See  Section  2.3  for  discussion.) 

4.  What  is  an  appropriate  initiation  lifetime  criterion? 

We  have  decided  on  equating  initiation  lifetime  with  the  cycles  to  form  a 
grain  sized  crack.  Our  reasoning  is  discussed  in  Section  2.3. 

5.  Finally,  how  can  one  simplify  the  complexities  of  nature  sufficiently  in  a 
model  to  allow  lifetime  to  be  predicted  under  spectrum  loading? 

At  first,  we  thought  that  it  might  suffice  to  model  what  happens  only  in 
the  largest  grains.  But  now  this  is  clearly  inadequate.  We  have 
therefore  designed  an  initiation  lifetime  model  that  accounts  for  grain 
size  effects.  Its  format  is  such  that  important  stochastic  elements  and 
perhaps  even  a  probabilistic  formulation  of  the  model  could  be  made  in 
the  future.  (See  Section  4  for  details.) 
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2.1  Mechanical  Property  Evaluation 

4  5 

In  two  previous  papers,  ’  we  used  models  of  the  deformation  of  a  soft 
ellipsoid  in  an  elastic  matrix  to  calculate  the  flow  stress  and  strain  hardening 
characteristics  of  surface  grains  in  A1  2219-T851.  We  found  that  fatigue  of  A1  2219-T851 
at  stresses  below  the  bulk  cyclic  yield  strength  caused  a  reduction  in  the  local  0.2% 
offset  yield  strength  of  the  interior  of  the  grains.  Grains  3-8  times  larger  than  the  60  ym 
mean  size  were  the  most  affected  and  are  the  eventual  sites  for  crack  initiation. 
Ultimately,  the  yield  strength  of  these  grains  falls  to  half  the  bulk  cyclic  value. 
However,  the  local  strain  amplitude  remains  small  because  the  deformation  is 
constrained  by  the  surrounding  nearly  elastic  grains. 

To  provide  more  experimental  evidence  for  this  fatigued  induced  reduction  in 
the  local  flow  stress,  we  used  a  load  reduction  sequence  to  minimize  the  influence  of 
local  residual  stresses.  This  technique  was  used  to  characterize  the  evolution  of  the  flow 
stress  in  a  300  pm  grain  in  A1  2219-T851.  The  results  were  published  in  Fatigue  87  (3rd 
Int.  Conf.  on  Fatigue  and  Fatigue  Thresholds)  (Ref.  6)  and  are  included  as  Appendix  A.  In 
summary,  two  flow  stresses  were  identified.  Microplastic  flow  takes  place  in  all  grains 
at  fully  reversed  stress  amplitudes  above  70  MPa,  although  the  plastic  modulus  is  large 
and  only  small  amounts  of  plastic  deformation  occur.  With  continued  fatigue,  an  upper 
flow  stress  develops  at  -  200  MPa  in  grains  substantially  larger  than  the  mean  grain 
size.  In  these  grains,  plastic  flow  is  easier  and  considerable  plastic  strain  occurs.  This 
suggests  the  existence  of  matrix  slip  (which  occurs  above  the  lower  flow  stress)  pene¬ 
trated  by  slip  bands  (which  occurs  above  the  upper  flow  stress).  The  deformation 
behavior  of  the  grain  might  be  explained  then  by  a  composite  model  incorporating  two 
deformation  regimes. 
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2.2  Yield  Locus  Motion  During  Spectrum  Loading 

The  behavior  of  the  two  stages  of  flow  identified  in  A1  2219-T851  was  investi¬ 
gated  in  more  detail  in  order  to  accurately  describe  the  deformation  pattern  by  a 
composite  model.  Various  load  sequences  were  used  since  prior  to  this,  all  our  studies 
had  been  done  using  only  fully  reversed  loading.  Two  limiting  cases  of  the  flow  surface 
locus  were  examined  by  studying  the  response  to  several  load  sequences.  These  were  the 
stationary  locus,  wherein  the  flow  stress  is  independent  of  prior  load,  and  the  kinematic 
surface,  in  which  the  flow  stress  moves  freely  with  the  maximum  stress  excursion  on  the 
previous  load  reversal.  This  effort  was  reported  at  the  Symposium  on  Effects  of  Load 
and  Thermal  Histories  on  Mechanical  Behavior  of  Metals^  and  is  included  as  Appendix  B. 

In  summary,  flow  above  the  lower  flow  stress  was  found  to  strain  harden 
quickly  and  is  best  described  by  kinematic  movement  of  the  flow  surface.  This  suggests 
there  is  considerable  reaction  to  short  range  back  stresses  in  the  grain,  and  deformation 
is  extremely  dependent  on  the  prior  load  sequence.  Deformation  above  the  upper  flow 
stress  shows  strong  strain  softening  with  fatigue  and  is  best  described  by  a  stationary 
yield  locus.  This  is  consistent  with  long  range  flow  constrained  by  grain  boundaries,  such 
as  in  coarse  slip  bands.  A  physical  representation  of  this  behavior  is  that  of  a  matrix 
containing  slip  bands.  A  uniform  deformation  model  incorporating  both  of  these  kinds  of 
behavior  was  developed  which  provided  an  adequate  description  of  deformation  in  the 
isolated  grains  responsible  for  crack  initiation. 
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2.3  Observations  Bearing  on  Slip  Band  Formation,  Hardening  and  the  Crack 

Initiation  Criterion 

Several  years  ago  we  showed,  for  A1  2219-T851  in  fully  reversed  loading,  that 

9 

low  amplitude  cycles  hardened  the  surface.  This  reduced  the  subsequent  initiation  of 
cracks  once  the  cyclic  stress  was  raised.  The  model  we  proposed  at  the  time  related  a 
hardening  rate  to  the  local  plastic  strain  range.  In  light  of  our  most  recent  findings, 
however,  the  resulting  expression  for  evolution  in  the  plastic  strain  range  cannot  be 
correct,  although  the  disparity  with  reality  might  be  difficult  to  notice  if  the  load 
spectra  is  not  complex.  The  old  model  in  effect  assumed  that  the  grain  hardening  was 
caused  by  cross-slip  which  interdicted  the  active  bands.  With  sufficient  fatigue,  the 
model  equations  predicted  that  a  grain  would  fully  harden  at  emy  cyclic  stress.  New 
measurements  of  strains  in  individual  grains  show  this  picture  is  inaccurate.  What 
actually  happens  is  that  raising  the  stress  after  hardening  causes  a  new  round  of 
softening,  although  producing  plastic  strain  amplitudes  lower  than  would  be  expected 
(Fig.  2.3-1)  if  one  had  begun  fatigue  with  the  higher  cyclic  stress.  Presumably, 
therefore,  it  is  only  the  slip  bands  formed  at  the  low  amplitude  which  harden. 
Apparently,  there  is  a  distribution  in  the  local  matrix  strain  range  required  for  band 
nucleation,  and  a  critical  strain  range  must  be  exceeded  to  produce  an  active  band  at 
each  potential  band  nucleation  site. 

How  can  we  relate  the  local  strain  range  in  the  matrix  and  the  strain  range 
required  for  development  of  a  band  at  each  potential  nucleation  site  in  a  surface  layer  to 
the  number  of  cycles  required  to  form  a  band?  Clearly,  there  is  insufficient  information 
about  the  microscopic  mechanisms  to  do  this  exactly,  so  we  make  the  following  rationali¬ 
zation.  For  A1  2219-T851,  the  fatigue  limit  in  fully  reversed  loading  is  ■  ±50%  of  the 
bulk  cyclic  yield  strength.  We  assume  that  the  plastic  matrix  strains  corresponding  to 
this  stress  range  are  just  of  sufficient  magnitude  to  produce  slip  bands.  For  applied 
stresses  greater  than  ±50%  of  yield,  the  number  of  cycles  (N)  to  produce  bands  is 
obtained  from  plots  of  the  local  residual  strain  vs  cycles,  such  as  Fig.  2.3-1.  The  cycles 
to  produce  a  band  at  a  site  requiring  the  maximum  applied  cyclic  strain  is  assumed  to 
correspond  to  the  cycles  to  reach  the  maximum  plastic  strain  amplitude.  (Because  each 
grain  cyclically  hardens,  this  procedure  probably  underestimates  N,  i.e.,  it  is 
conservative.)  The  cycles  to  produce  bands  at  sites  capable  of  forming  bands  at  50%  of 
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Fig.  2.3-1  Residual  strain  measured  in  a  300  um  grain  after  tensile  loading. 

Specimen  first  cycled  until  hardened  at  ±0.8  Oy,  then  subjected  to 
±0.9  Oy.  [oy  is  bulk  yield  strength.]  ^ 


yield,  when  the  applied  stress  is  greater  than  50%,  is  read  as  the  cycles  where  the  onset 
of  a  positive  residual  strain  begins.  From  Fig.  2.3-1,  N  s  150  cycles  at  ±0.8  Oy. 

We  propose  an  accumulative  damage  model  to  predict  the  cycles  to  form  each 
slip  band.  This  is  calibrated  for  constant  stress  amplitude  (and  quasi  constant  local 
matrix  strain  amplitude)  loading. 

Imagine  that  a  grain  experiences  a  local  plastic  strain  range  Ae^  in  the 

matrix.  The  total  measureable  strain  increases  with  fatigue  as  slip  bands  form,  until  the 

band  initiation  process  stops.  At  this  point,  the  range  Ael"  has  activated  all  bands 

I  ^ 

requiring  an  activation  strain  range,  Ae^,  which  is  less  than  the  plastic  strain  in  the 
matrix,  Ae™  .  No  further  bands  can  be  activated  unless  A£p  is  increased. 
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However,  before  this  point  is  reached,  all  sites  requiring  activation  strains 
acp  <  ACp  have  produced  bands.  A  simple  construction  described  below  allows  us  to 
estimate  the  number  of  constant  amplitude  cycles  this  takes.  Later  on,  as  we  develop 
the  spectrum  loading  model,  the  reciprocal  of  these  cycles  (£q.  2.3-2)  will  be  the  damage 
rate  for  the  local  plastic  strain  range  Ac^  . 

We  can  use  an  expression  for  the  local  fully  reversed  plastic  strain  range  in  a 
kinematic  domain  (the  matrix  surface  layer  denoted  by  m)  embedded  in  the  grain  to 
relate  the  strain  range  to  Og,  the  local  stress  excursion  in  the  grain. 


(2.3-1) 


Oj  has  already  been  defined  as  the  lower  flow  stress,  is  the  plastic  modulus  in  the 
surface  layer,  a  is  a  shape  factor  constant  and  Eg  is  Young's  modulus.  This  allows  us  to 
make  the  construction  in  Fig.  2.3-2.  The  symbols  represent  data  points  from  experiments 
where  the  cycles  required  for  the  local  stress  to  exceed  and  were  obtained  at  four 
different  applied  stress  levels.  Once  the  local  stress  exceeds  o^,  the  plastic  strain  in  the 
matrix,  A£p  ,  exceeds  that  for  activation  of  slip  bands.  The  solid  line  in  Fig.  2.3-2  is  fit 
through  the  second  of  each  pair  of  data  points  taken  at  different  applied  stress  ampli¬ 
tudes.  This  line  relates  the  critical  plastic  strain  necessary  to  activate  all  slip  bands  to 
cycles.  Some  slip  bands  activate  at  a  lower  strain  down  to  that  reached  when  is 
exceeded,  as  denoted  by  the  lower  data  point  of  each  pair.  The  dashed  line  connecting 
the  pair  of  experimental  data  points  to  the  origin  can  then  be  used  to  calculate  the 
cycles  required  to  form  a  band  at  a  site  requiring  a  lower  activation  plastic  strain  than 
the  applied.  Call  Acp  the  local  range  for  band  activation  and  the  cycles  to  band 
formation  seen  in  Fig.  2.3-2  has  the  form 


(2.3-2) 
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Yj  and  Y2  material  parameters  determinable  by  fitting  Eq.  (2.3-2)  to  the  data  in 
Fig.  2.3-2.  We  have  noted  earlier  in  this  report  that  deformation  below  the  upper 
flow  stress,  is  kinematic  and  is  therefore  in  the  matrbc  penetrated  by  bands.  As  such,  the 
local  deformation  is  short  ranged  and  can  "feel"  the  presence  of  grain  boundaries  only 
through  a  reaction  stress  to  the  net  plastic  strain.  Thus,  n  in  Eq.  (2.3-2)  is  independent  of 
grain  size. 


CONSTANT  AMPLITUDE  1/2  CYCLES 


Fig.  2.3-2  Relationship  of  cycles  to  band  initiation  to  Acp  and  Aep  for  fully 

reversed  loading.  The  dashed  lines  represent  families  of  critical  plastic 
strain  to  band  activation  at  constant  plastic  strain  in  the  matrix  surface 
layer. 


However,  this  does  not  mean  that  the  strains  measured  in  each  grain  will  be 
independent  of  grain  size.  Observations  made  under  a  prior  contract  with  NADC  have 
led  to  the  hypothesis  that,  in  A1  2219-T851,  slip  is  initiated  in  a  thin  layer  near  the 
surface.  The  number  of  band  nucleation  sites  in  each  grain  will,  therefore,  be 
proportional  to  the  grain  width  (0)  and  the  band  displacements  per  unit  length  will  be 
proportional  to  0,  giving  «  0^  . 
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Such  a  trend  has  been  verified  for  Al  2219-T851  using  data  acquired  several 
years  ago  (Fig.  2.3-3).  Additionally,  another  set  of  residual  strains  in  150  and  300  pm 
grains  have  been  determined  for  Al  2219-T851.  These  were  taken  after  several  hundred 
cycles  and  also  confirm  trend  with  grain  size  shown  in  Fig.  2.3-3.  So,  the  interiors  of  all 
the  surface  grains  apparently  soften  at  the  same  rate  in  terms  of  bands  produced  per  site 
available,  but  the  grain  size  affects  the  magnitude  of  the  strains  which  result. 


0  1  00  200  300  400  600 

SLIP  DISTANCE  Urn) 


Fig.  2.3-3  Effect  of  grain  size  on  measured  residual  strain.  Data  from  Ref.  9  for 

±270  MPa  cycles  in  Al  2219-T851.  Residual  strains  measured  after  a 
tensile  load.  Curves  are  proportional  to  and  fit  the  data  for  grains 
small  enough  to  be  essentially  unhardened. 


The  other  important  effect  of  grain  size  appears  in  the  grain  hardening  rate. 
We  think  this  is  because  the  displacements  experienced  by  individual  bands  increase  with 
grain  size.  Large  grains  harden  much  faster  than  small  grains.  Also,  there  is  a 
pronounced  effect  of  applied  strain  range  on  the  rate  of  hardening  of  individual  grains 
(Fig.  2.3-4).  The  hardening  rate  decreases  markedly  as  the  cyclic  stress  is  lowered. 
Because  the  displacement  of  individual  bands  should  be  proportional  to  (Aa^  -  2O|^)0  in 
fully  reversed  loading,  there  should  be  a  correlation  of  this  quantity  to  the  cycles  for 
band  hardening.  We  use  the  number  of  constant  amplitude  cycles  it  takes  the  locally 
measured  residual  strain  to  fall  to  ■  1/e  of  its  maximum  value  to  obtain  an  empirical 
relationship  of  the  cycles  to  harden  a  band  as  a  function  of  the  grain  size  and  load 
(Fig.  2.3-5).  In  fully  reversed  loading. 
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Fig.  2.3-4  Effect  of  loading  range  on  hardening  of  residual  strains  measured  in  a 
300  ym  grain.  Compare  to  Fig.  2.3-1  where  hardening  at  constant 
amplitude  is  almost  complete  by  600  cycles. 


•C4t191 


Fig.  2.3-5 


Empirical  relationship  of  estimated  cycles  to  fall  in  «  i/:' 
maximum  local  residual  strain  and  the  quantity  (ac^  -  2o 
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N  =  10 


3 


-  20^)0) 


(2.3-3) 


exp>resses  the  relationship  in  Fig.  2-3.5.  The  t's  are  constants  determinable  for  A1  2219- 
T851  from  Fig.  2.3-5.  Once  a  band  forms,  it  presumably  forgets  the  Acp  of  its  initiation 
site,  and  only  the  strains  to  which  it  is  exposed  after  it  forms  and  their  duration  control 
its  hardening  state. 

We  now  discuss  our  rationale  for  defining  crack  initiation  as  the  cycles  to 
obtain  a  crack  of  one  grain  depth.  Experiments  we  did  several  years  ago  showed  a 
strong  correlation  between  the  peak  local  plastic  strain  in  a  grain  and  constituent 
particle  fracture.  However,  those  grains  which  deformed  extensively  also  accumulated 
large  cycle-by-cycle  damage  levels,  casting  doubt  as  to  the  true  initiation  mechanism, 
i.e.,  peak  strain  or  accumulated  damage.  Furthermore,  for  this  study,  there  was  the  need 
of  a  deformation  model  of  initiation  with  sufficient  practical  engineering  significance  to 
formulate  a  lifetime  model.  Since  single  extreme  overloads  can  fracture  constituent 
particles  in  Al  2219-T851,  but  not  propagate  the  crack  into  the  grain,  a  peak  strain 
criterion  describing  this  behavior  seemed  of  no  practical  importance.  However,  early 
crack  growth  to  the  first  grain  boundary  comprises  a  phenomena  common  to  many  alloys, 
including  Al  2219-T851,  which  is  also  a  type  of  "initiation"  and  can  be  a  primary  factor  in 
controlling  lifetime.  Our  modeling  of  local  plasticity  provides  for  the  first  time  the 
opportunity  to  predict  this  "initiation  liie". 

For  15  years,  it  has  been  recognized  that  very  short  surface  cracks  propagate 

11  12 

faster  than  long  cracks  at  the  same  cyclic  stress  intensity.  After  corrections  have 

been  made  for  closure  stress  and  its  dependence  on  crack  length,  and  for  crack  growth 
statistics  affected  by  local  microstructure,  one  finds  that  the  very  shortest  cracks 
(perhaps  smaller  than  25  ym)  still  grow  much  too  fast  to  be  understood  by  linear  elastic 
fracture  mechanics.  Lankford  et  al,^^  Smith*^  and  Miller^^  have  attributed  this 
accelerated  growth  to  an  enhanced  crack  tip  plasticity,  perhaps  predictable  by  a  second 
order  term  in  the  expression  for  stress  intensity.  However,  for  very  short  cracks,  these 
stress  intensity  corrections  still  grossly  underestimate  the  actual  growth  rates.  This,  we 
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believe,  is  because  these  models  presume  that  the  surface  is  otherwise  elastic,  while  it  is 
clearly  not  in  many  materials.  In  prior  work,^^  we  found  that  by  using  the  elastic-plastic 
crack  growth  expression  by  Tompkins,^  ^  we  obtained  from  known  plastic  strains  in  a 
grain  at  a  crack  tip  an  estimate  of  the  true  crack  tip  opening  displacement  (CTOD)  for 
cracks  50  ym  and  larger.  It  was  typically  a  factor  of  two  larger  than  anticipated  for  an 
elastic  medium.  It  was  the  long-range  plastic  strains  induced  not  by  the  crack  itself,  but 
by  the  independent  deformation  of  the  grain  surrounding  the  crack  which  increased  the 
CTOD  and  accelerated  growth. 

Further  analysis  now  shows  that  the  relatively  small  correction  to  CTOD 
provided  by  the  Tompkins  expression  may  produce  several  orders  of  magnitude  increase 
in  growth  rate  to  a  microcrack  smaller  than  20  ym,  the  typical  size  of  the  largest 
intermetallic  in  A1  2219-T851.  The  key  in  predicting  growth  rate  is  to  accurately  know 
the  local  plastic  strain  range  cycle-by-cycle.  In  all  likelihood,  as  long  as  the  entire  crack 
is  embedded  in  the  plastic  strain  field  of  a  surface  grain,  it  will  feel  the  accelerating 
effect  of  this  plasticity.  Once  a  crack  begins  to  pass  below  the  first  subsurface  grain 
boundary,  the  growth  rate  should  start  to  converge  towards  an  LEFM  prediction.  This 

t  0 

association  was  noted  by  Lankford  some  years  ago,  but  the  reason  why  it  existed  was 
mystifying.  Under  this  contract,  we  have  measured  strains  on  the  surface  just  in  front  of 
a  growing  15  ym  crack  and  confirm  that  the  grain  must  be  softer  for  there  to  be 
substantial  enhancement  propagation  in  A1  2219-T851. 

This  simple  definition  of  initiation,  as  the  number  of  cycles  required  to  initiate 
and  propagate  a  microcrack  to  the  first  subsurface  grain  boundary,  has  satisfying  engi¬ 
neering  applications.  The  arrival  of  the  crack  front  at  the  subsurface  grain  boundary 
corresponds  to  the  point  in  lifetime  beyond  which  LEFM  suffices  to  give  a  crude  estimate 
of  the  life  remaining  (although  corrections  for  effects  of  the  microstructure  on  the  small 
crack  growth  would  still  be  advisable).  This  is  not  to  be  considered  a  universal  definition 
of  "initiation,"  but  it  should  apply  to  many  high-strength  aluminum  alloys. 
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3.0  APPARENT  "INITIATION"  LIFETIME  SCENARIO 

Irwiividual  grains  develop  slip  bands,  whose  presence  strain  softens  the  local 
surface,  accelerating  crack  initiation.  In  A1  2219-T851,  the  bands  are  thought  to 
nucleate  at  sites  scattered  within  a  thin  («  2  wm  deep)  surface  layer,^  as  illustrated  in 
Fig.  3.0-1.  For  a  fully  reversed  local  plastic  strain  in  the  layer,  damage  leading  to  band 
nucleation  accumulates  faster  the  larger  the  strain  &ep  and  the  smaller  the  activation 
strain  range  (Acp  )  for  each  site.  &cp  must  be  greater  than  or  equal  to  Acp  for  a  band  to 
ultimately  form  by  fatigue.  A  phenomological  model  for  the  rate  of  band  initiation 
vs  Acp  and  Atp  has  been  proposed  in  Section  2.3  for  A1  2219-T851,  based  upon 
microscopic  strain  observations.  Since  deformation  of  the  surface  layer  is  kinematic, 
there  can  be  no  direct  effect  of  grain  size  on  Acp  .  However,  the  layer  size  affects  the 
number  of  slip  bands  that  can  form,  since  larger  grains  will  have  more  potential 
nucleation  sites/length.  Also,  the  coupling  between  grain  size,  band  density  and  elastic 
modulus  in  the  grain  body  will  affect  the  local  stress  range  which  controls  Aep  . 


•C41129 


Fig.  3.0-1  Schematic  cross-section  of  a  grain  showing  bands  that  have  nucleated  at 
sites  in  a  thin  surface  layer. 


18 

C8577D/bje 


NADC-89044-60 


Rockwell  International 

Science  Center 


SC5470.fr 


To  find  Atp  and  the  plastic  strain  range  in  the  grain  body  Ae^  ,  it  is  necessary 
to  make  a  parallel  calculation  of  the  deformations  in  the  grain  body  vs  band  density,  as 
well  as  in  the  surface  layer.  For  economy  of  numerical  analysis,  the  load  spectrum  will 
be  subdivided  into  repetitive  blocks. 


3.1  Efface  Layer  -  Nucleated  Band  Density 

An  incremental  change  in  damage  at  a  potential  band  nucleation  site  over  a 
loading  block  is  defined  to  be  Al(Epg,  Acp).  [Naturally,  Al  also  depends  upon  the  load 
spectra  comprising  the  block.]  Ep^  is  the  plastic  modulus  of  the  entire  grain.  Summed 
over  repeated  blocks,  the  requirement  for  initiation  of  a  band  at  a  site  of  initiation 
range  Acp  is 

I(AeJ)  =  1  .  (3.1-1) 

After  m  blocks,  there  exists  a  maximum  activation  strain  of  the  sites  which  have 

produced  slip  bands.  We  call  this  quantity  and  it  is  found  by  solving  Eq.  (3.1-1). 

If  no  strain  hardening  has  occurred  during  m  blocks,  the  linear  density  of 
active  slip  bands  p  in  a  grain  of  width  D  would  simply  be 

0  =  0  J  ^  Ti(Aep)dAep  .  (3.1-2) 

0 

ri(Aep)  is  the  number  of  band  nucleation  sites/area  of  each  type. 

3.2  Grain  Body 

Under  stress,  each  active  band  will  displace  an  amount  proportional  to  its 
width.  The  apparent  plastic  modulus  of  the  grain  becomes 

r  ^0 

n(..‘)  d..‘ 
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9 

b 


(3.2-lb) 


where  Eq  and  E^'  are  material  parameters. 

E--  is  an  effective  modulus,  i.e.,  an  average  of  the  composite  caused  by 

•  O 

plastic  flow  of  the  matrix  and  slip  of  any  active  bands  present.  When  there  are  no  bands, 
Ejj  •  and  Epg  =  E^',  which  must  therefore  be  the  plastic  modulus  of  the  matrix,  i.e., 
Eq'  "  ^pm‘  Parameter  E^  can  be  determined  by  calibration  once  Eq.  (3.2-la)  is 
calculable.  This  requires  a  determination  of  Epg(n)  as  described  in  Sect.  2.2  (defined  as 
Ep(i)  in  the  paper)  for  comparison  to  Eq.  3.2-lb  through  the  application  of  Eqs.  3.2-la  and 
2.3-2.  Since  Eq  is  a  fundamental  material  parameter,  this  analysis  needs  to  be  done  just 
once  for  each  material  modeled. 

When  a  spectrum  load  is  applied,  the  local  stress  in  the  grain's  main  body 
Og  differs  from  the  applied  by  a  residual  stress  which  is  proportional  to  the  local 
plastic  strain  in  the  grain  (e^ ).  The  stress-strain  behavior  is  calculated  by  a  "stationary" 
flow  criterion,  which  treats  a  single  flow  stress  of  the  active  bands  as  invariant. 

Softening  in  this  model  is  the  result  of  a  progressive  reduction  in  Ep^  as 
increases. 

Strain  hardening  is  handled  by  considering  the  deformation  of  each  slip  band. 
We  assume  in  accord  with  the  experiment  that  each  new  band  will  slip  on  exceeding 
until  it  has  accumulated  a  hardening  level  between  m’(nucleation)  and  blocks 

-  I-  (3-2-2> 

The  cycles  to  reach  this  critical  ^  will  depend  upon  grain  size  because  the  cycle- 
by-cycle  strain  in  each  individual  band  increases  with  grain  size  for  a  constant  Epg. 

We  define  an  incremental  hardening  aH (Epg,  D)  for  one  loading  block, 
determined  by  microscopic  observations.  (Epg  and  p  depends  upon  0.  )  aH  is  summed 
block-by-block  to  determine  aHj^i^pp  and  thus  how  long  it  takes  a  newly  formed  band  to 
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cease  slipping.  Once  a  band  forms,  it  apparently  loses  all  memory  of  its  origin  other  than 
•Vhen"  it  formed.  All  bands  nucleated  at  m*  blocks  will  therefore,  in  this  model,  cease  to 
slip  at  the  same  time  in  the  future.  So,  the  density  of  active  slip  bands  (p)  can  be 
expressed  by  a  simple  integral  of  n* 


where 


I 


max 


.4ep 

=  0  J  j  n(Aep)  dAcp 

,  ^min  ^  ^ 


.  ^mln  _  ^ 

AEp  =  Acp  (m  ):  and 


m'  is  obtained  by  solving 


(3.2-3) 


(3.2-4) 


H(m)  -  H(m')  =  1. 


(3.2-5) 


Wherefore 


I 

,2  r 


max 


0  J  j  n(Aep)  dAtp 

min 


■'  -cl 


Ac, 


-  ^o'^b 


(3.2-6a) 


(3.2-6b) 


Consequently,  grain  hardening  with  fatigue  in  this  model  is  tied  to  a  progressive  increase 

in  Acp'"^^  the  minimum  plastic  strain  range  capable  of  initiating  a  slip  band,  and  the 
resulting  increase  in  Ep^. 


3.3  Crack  "Initiation" 

First,  we  must  calculate  Ep^  and  the  associated  and  cycle-by-cycle. 
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These  stresses  and  strains  control  initiation.  As  a  practical  matter,  the  smallest  cracks 
that  are  routinely  found  in  laboratory  specimens  in  structural  aluminums  are  typically 
one-grain  thickness  in  depth.  Since  local  plasticity  in  a  surface  grain  greatly  accelerates 
growth,  a  grain-sized  crack  corresponds  to  the  smallest  crack  size  for  which  simple 
microstructual-based  corrections  to  linear  elastic  fracture  mechanics  provide  an 
adequate  expression  for  growth  rate.  Also,  the  average  growth  rate  is  a  minimum  at 
one-grain  depth  (Fig.  3.3-1). 


Fig.  3.3-1  Typical  dependence  of  average  short  crack  growth  rate  on  aK. 

Plastic  strains  can  accelerate  growth  in  the  first  grain  by  a  factor  of  100.  So, 
a  definition  of  "initiation”  as  corresponding  to  one-grain  depth  has  fundamental  physical 
and  engineering  significance.  In  the  following,  we  show  how  our  current  insight  on  local 
plasticity  can  be  incorporated  into  an  elastic-plastic  growth  equation  of  Tompkins (see 
Section  4.1.7)  to  predict  "initiation”  under  spectrum  loads.  We  distinguish  Mode  I  crack 
growth  initiating  from  intermetallic  particles  and  mixed  mode  growth  from  slip  bands,  so 
that  the  model  should  be  suitable  for  most  high  strength  aluminum  alloys  having  these 
crack  initiation  modes. 
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4.0  OVERALL  MATHEMATICAL  STRUCTURE 

The  cycles  to  initiation  (mn)|  [U  of  blocks  times  //  cycles  in  block]  will  depend 
on  the  local  grain  size.  The  grain  size  leading  to  earliest  initiation  will  vary  with  the 
load  spectrum.  Brute  force  numerical  simulations  of  initiation  must  therefore  be 
replaced  with  a  mathematics  of  the  initiation  process  suitable  for  solution  by  numerical 
techniques.  The  numerical  coding  of  the  response  functions  needed  to  exercise  this 
model  is  left  to  the  reader.  In  the  following,  the  model  is  stated  in  detail  only  for 
initiation  with  mode  I  growth  from  constituent  particles. 

Because  the  rate  at  which  the  local  mechanical  properties  change  depends  on 
grain  size,  we  cannot  a  priori  tell  which  local  grain  size  in  an  alloy  will  lead  to  earliest 
initiation.  In  fact,  it  is  certain  that  this  critical  size  will  depend  upon  load  sequence. 
The  trend  will  be  for  large  grains  which  can  sustain  layer  plastic  strains  to  advance  a 
microcrack  most  rapidly  from  a  fractured  intermetallic.  However,  large  grains  are  also 
more  subject  to  strain  hardening,  so  they  may  only  actively  deform  for  a  short  time, 
leaving  the  job  of  producing  the  terminal  crack  to  a  grain  of  more  modest  size.  In  the 
following,  we  define  a  set  of  damage  functions  which  depend  on  grain  size.  The 
numerical  concept  is  to  determine  these  functions  for  block -by -block  loading,  so  that  the 
grain  size  dependence  of  the  incremental  crack  growth  is  automatically  contained  in 
each  intermediate  result.  The  calculation  can  be  carried  out  until  the  crack  depth  within 
the  range  in  grain  sizes'  reaches  a  critical  size.  The  calculation  does  not  follow  the 
evolution  of  damage  of  single  grains,  but  rather  describes  the  evolution  of  a  damage 
state  space  in  which  grain  size  is  a  variable. 

Four  levels  of  analysis  are  identified.  At  level  1  are  three  individual  loading 
block  response  functions  in  which  the  response  of  any  surface  grain  is  embedded,  at  any 
time  during  fatigue,  to  one  standard  loading  block.  The  parameters  represented  will 
allow  the  progressive  evolution  of  the  local  mechanical  properties  to  be  predicted,  from 
which  the  incremental  change  in  crack  length  of  a  subgrain-sized  crack  can  be  estimated 
block-by-block. 
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The  functions  are; 


‘KEp,.  A.J)  i 

incremental  damage  for  band 
nucleation 

AHIEpg.  0)  i 

incremental  damage  for  band 
hardening 

Aa(Ep^.  a)  5 

irKremental  growth  of  a 
subgrain  sized  crack  of  depth 

The  new  variable  here  is  (a),  the  crack  depth.  The  way  these  functions  serve  the 
initiation  lifetime  calculation  is  illustrated  schematically  in  Fig.  4.0-1. 

Once  these  three  functions  have  been  determined,  a  simulation  of  the 
initiation  process  can  proceed.  The  functions  will  be  used  as  "look  up"  tables  for  values 
needed  in  this  sequential  analysis. 


•C411»3 


Fig.  4.0-1  Block-by-block  calculation  of  crack  length. 
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At  level  2  are  functions  which  must  be  determined  and  stored  only  to  the  next 
loading  block.  These  are 


UD, 


the  cumulative  damage  for  initiation  of  a 
band  at  a  Acp  activation  site  in  a  grain  of 
size  D 


a(D); 


the  crack  depth  vs  D 


max /n \ , 
aeq  (0); 


the  maximum  &Ep  which  has  caused  band 
nucleation 


aH(D)  ;  the  incremental  hardening  vs  grain  size 

Two  level  3  functions  retain  a  running  memory  from  the  first  block  of  the 
cumulative  damage  governing  band  density.  These  are  necessary  to  calculate  the 
mechanical  property  state  for  each  succeeding  block. 


tt 


^max 

P 


(D.m)  ; 


is  a  running  value  of 


Ae 


max 


(D) 


H(D,m) ;  is  the  summed  hardening  coefficient 

The  final  level  4  contains 

Aep"'^^(0)  ;  the  minimum  Acp  for  which  bands  have 

ceased  to  slip 

Ep^;  the  plastic  modulus  of  the  grain's  main  body 

p(D) ;  the  active  slip  band  density 
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Some  reasonable  mixture  of  microscopic  data  on  local  strains  and  fatigue  test 
data  will  be  needed  to  calibrate  this  model.  The  following  is  just  one  of  several  options. 

4.1  The  Loading  Block  Response  Functions 

The  spectrum  loading  block  must  first  be  expressed  in  terms  of  sequential 

minimum  and  maximum  applied  stresses  and  ®  selected  Ep^,  these  are 

used  to  find  the  local  stress  and  plastic  strain  in  the  grain  body,  is  then  used 

as  an  effective  applied  stress  to  find  Acp  in  the  surface  layer  embedded  within  the 

19  a 

grain.  With  appropriate  rain  flow  sequence  counting  acq  and  Ao  give  Aa,  Ao  is  used 

fn  ^  ^ 

to  find  aH,  and  Acp  is  used  along  with  ^pg(D)  4l.  See  Fig.  4.1-1  for  the  logic  of 

this  procedure. 


Fig.  4.1-1  Calculating  the  loading  block  response  functions. 

Note;  The  calculation  described  is  of  an  average  initiation  lifetime. 
Variability  in  constraint  (a^  and  for  the  surface  layer  and  grain,  respectively)  and  in 
o^j  from  grain-to-grain  could  be  included  to  treat  statistical  scatter. 
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Rain  flow  counting  techniques  (for  which  "real-time"  algorithms  are 
available)^®  will  be  used  to  express  the  local  strains  and  stresses  found  by  our 
deformation  equations  in  terms  of  equivalent  closed  cycles.  Five  cumulative  probability 
distributions  describing  these  quantities  are  needed  (see  Fig.  4.1-2  for  an  example). 


■caiist 


P(o.?(Ep,I)i 

for  the  plastic  strain  range  in  the  grain  body 

P(‘«glEpgl)i 

for  the  stress  range  in  the  grain  body 

P('‘^"glEpgl)i 

for  the  positive  stress  range  in  the  grain 
body 

P(‘^?lEpgl)i 

for  the  positive  plastic  strain  range  in 
the  grain  body 

P(“>?|Epgl)i 

for  the  plastic  strain  range  in  the  surface 
layer 
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denotes  the  tensile  portion  of  the  cycle  only,  since  crack  growth  only  occurs  on 
positive  loading.  The  effects  of  load  sequence,  ntean  stress  and  stress  amplitude  are 
embedded  in  these  distributions  through  their  affect  on  the  local  stress-strain  behavior. 


4.1.1  Rain  Flow  Couiting  and  the  P  Distributions 

A  typical  output  from  the  cycle-by-cycle  strain  calculations  described  in 
Sections  4.1.2  and  4.1.3  will  be  as  in  Fig.  4.1-3a.  Rain  flow  analysis  will  be  needed  to 
reduce  this  spectra  to  equivalent  closed  cycles.  The  scheme  is  illustrated  schematically 
in  Fig.  4.1-3a  and  the  result  is  tabulated  in  Fig.  4.1-3b  giving  a  P  distribution  for  one 
value  of  Epg.  In  the  following  subsections,  we  discuss  how  the  raw  spectra  on  which  this 
procedure  will  operate  are  obtained,  and  then  how  the  block  response  functions  are 
calculated. 

4.1.2  The  Statiortary  Grain  Body  Deformation  Model  Uog,  and  &e^] 

Assume  o*,  and  of are  the  applied  maximum  and  minimum  stresses.  There 
max  mm 

are  five  cases: 

Case  I:  Loading  from  to  ,  initial  Og  exceeds  o^j. 


I  O  -  0  .  -  o  +  o  1 

‘max  mm  u  q 

=  - F - 7 - C -  ■  On 

3  pg  g  e 


(4.1-1) 


In  -O  .  -  0  +  ol 

.  q  ‘  max  mm  u  q' 

"  E  +  a  E 

pg  g  e 


(4.1-2) 


where  is  the  bulk  elastic  modulus. 


Case  II;  Loading  from  o^^^^  to  with  initial  Og  >  0^^ 


Ao. 


^^max  ~  °min^^pq 
^pg  “g^e 


(4.1-3) 
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tC41132 

a) 


e 


TIME 


Rain  flow  analysis  procedures  of  local  plastic  strain  on  stress  spectrum, 
(a)  Raw  spectra  and  rain  flow  lines  (dashes);  and  (b)  cumulative  spectra 
distribution. 
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.4. 


r  &  ^1 

^°max  ~  °min‘ 

E„_  +  o  E 
pg  g  e 


(4.1-4) 


Case  III:  Always  elastic. 


AOg  =  Ao 


(4.1-5) 


Acp  =  0 


(4.1-6) 


than  -o^J. 


Case  IV:  Unloading  from  to  0^^^  with  initial  >  -0^^  ,  final  load  less 


( +0^  +0  +  0  ] 

‘min  max  u  q  1  .  1 

AO  -  - F - r - ^  -  lo  +0 

3  pg  g  e 


(4.1-7) 


0  .  +0  +  0  +  0  1 

.  q  _  ‘  mm  max  u  q* 

Sg  *  »9‘e 

Case  V:  Unloading  from  0®, „  to  with  <  -0.. 

°  max  mm  g  u 


(4.1-8) 


I’mln  -  ‘’maxl^DQ 

9  “"TTTTT— “ 
»  pg  g  e 


(4.1-9) 


Ac 


r  a  a  1 

g  ^  *°m1n  ~  °  max* 

E  +  o  E_ 

pg  g  e 


4.1.3  Kinematic  Matrix  Layer  Deformation  Model  (Aep) 

Maximum  and  minimum  values  of  0^  taken  from  Sect.  4.1.2  are 

m  9  9  9 

used  to  find  Aep  .  There  are  three  cases.  Separate  values  for  the  upper  and 
lower  Oj^2  stress  of  the  layer  must  be  calculated  as  these  will  change  cycle-by-cycle 


(4.1-10) 
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because  the  flow  locus  is  kinematic.  The  local  stress  in  the  matrix  layer,  will  also 
change  from  the  local  stress  in  the  grain  body,  Og. 

Case  I:  Loading  from  to  Og”***  exceeding  the  prior  . 


,  max  min  ^ 

-  °q  “ 


ir 


E„_  +  a_E- 
pm  me 


(4.1-11) 


40^1= 


C  A  ^ 

pm  p 


(4.1-12) 


4a„  =  E„_Ae|!I'  +  o., 
m  pm  p  ti 


m 


(4.1-13) 


Case  II:  Always  elastic. 


40,1  =  40,2 


=  0 


(4.1-14) 


4£p  =  0 


(4.1-15) 


4o_  =  Ao- 
m  g 


(4.1-16) 


Case  111:  Unloading  from  Og"***  to  Og  ,  falling  below  the  prior  0,2* 


,  min  max  .  1 

o„  -  o_  -  0,0  * 
m  3 _ 3_ 


4eo  = 


E„+  a_E_ 
pm  me 


(4.1-17) 
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(4.1-18) 

(4.1-19) 


4.1.4  Slip  Band  Density,  n 

A  crude  estimate  of  n  can  be  made  from  observations  (Fig.  4.1-4)  of  the  peak 
strains  in  grains  vs  the  cyclic  stress  amplitude  for  fully  reversed  loading.  Our 
observation  are  consistent  with 

n  =  AeJ  .  (4.1-20) 

A  proportionality  constant  needed  in  this  expression  has  been  embedded  in  and  E^'  in 
the  equations  used  to  calculate  Ep^  and  p.  Values  at  larger  strains  are  taken  from 
specimens  first  hardened  at  smaller  amplitudes. 


(C41<]4 


Fig.  4.1-4  Estimate  of  n  made  from  peak  residual  strain  amplitudes. 
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4.1.5  Incremental  Hardening  [aH] 

Experiments  indicate  that  at  constant  local  cyclic  stress,  the  cycles  to  harden 
a  slip  band  can  be  expressed  by  Eq.  (2.3-3).  and  T4  are  material  parameters  found  by 
microscopic  deformation  (^servations.  For  spectrum  loading,  we  propose  a  linear 
hardening  rule,  and  obtain  from  Eq.  (2.3-3)  an  incremental  hardening  of 


10 


(Y4(40g  -  2o^)D} 


(4.1-21a) 


for  one  cycle. 

From  whence, 


ta  dP(Ao  )  1^4(60  -  2o  )] 

aH(E  ,0)  =  n  J  9  -5-^  10  ^10  ’  9  “  dio  .  (*.l-2lb) 

2  a.  ''  ^ 


An  alternate  approach,  perhaps  of  use  for  treating  microcracking  of  the  slip  bands 
themselves,  might  be  to  relate  AH  to  the  displacements  in  each  band.  These  will  be 
proportional  to  Ae^/p . 

4.1.6  Increment  in  Slip  Band  Nucleation  Damage  (a1) 

Apparently,  a  critical  strain  range  must  be  exceeded  in  the  surface  layer  for 
fatigue  to  activate  slip  at  a  potential  nucleation  site.  Furthermore,  these  critical  values 
of  Acp  vary  from  site  to  site.  Our  phenomological  model  which  relates  Acp  at  the  site  of 
activation  Acp  to  the  cycles  to  nucleate  a  band  at  constant  amplitude  is  given  in  Fig. 
2.3-2  and  has  been  calibrated  (yj,  y2)  by  microscopic  deformation  data  for  A1  2219-T851. 

We  assume  for  spectrum  loading,  a  linear  accumulation  of  damage,  c,  and 
deduce  for  one  cycle  that 
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(4.1-22) 


‘I(Epg.  »ej)  = 


(*.1-23) 


4.1.7  Irxrrement  in  Crack  Length  (Aa) 

First,  we  truncate  the  (jg  and  spectra  to  their  positive  values  before  the 
rain  flow  analysis  is  done.  Thereafter,  we  have 


P(A^Og)  and  P(a^9) 


containing  n^  and  Op  cycles,  where  n^  is  the  number  of  cycles  in  the  stress  spectra  and 
Op  is  the  number  of  cycles  in  the  plastic  strain  spectra,  respectively.  Because  some  of 
the  cycles  are  elastic. 


By  manipulating  a  model  for  crack  tip  opening  (aa)  by  Tompkins,  we  obtain  a  cycle- 
by-cycle  expression  for  the  elastic -plastic  incremental  extension  of  a  crack  of  length  a. 

Aa  *  A(A6)^n  (4.1-24) 


where 
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(4.1-25) 


from  adaption  of  Tomkins  equation.  Parameters  A  and  y  can  be  determined  by 
calibration  from  initiation  lifetime  data. 

Since  is  monotonically  related  to  we  can  perform  a 

renormalization  to  allow  Eq.  (4.1-25)  to  be  integrated  over  a  loading  block,  without 
relying  on  a  cycle-by-cycle  tracking  of  the  relationship  of  the  local  stresses  and  strains. 
The  complexity  to  be  cared  for  is  that  some  of  the  cycles  are  entirely  elastic. 

For  stresses  less  than  Og  defined  by  solving 


n.  -  n 


(4.1-26) 


there  is  no  plastic  strain. 

For  larger  stresses,  the  equivalent  A'e^  is  found  by  forcing  a  mapping  of  the 
stress  range  4'Og  -  AOg  to  the  strain  range  0  -  Ae^.  We  have 


^  iP(j.g)  -  P(«  'Ogjl  =  PlJ’.S), 


(4.1-27) 


and 


-  p'‘  (P(p\)  -  p(«^'p»j)i 


(4.1-28) 


Hence,  the  incremental  change  in  crack  length  equals  the  sum  of  the  elastic  and  plastic 
growth  terms. 
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J'o  dP(J.  )  a(l.l)^(4V)2’'  ^ 

*•  -  (n,  -  np)A  ;  5  -1—^  I - JE-^I  °g 

0  dA  Og  e  ^ 

A^o  dP(A^0  )  J  2  Y  (l.l)^A^o  T  'o  ^  T 

+  n  A  J  x  ®  — T  ^  l4(A  o  )  1  1 - Rp - ^  +  0.4»a  E^j  dA  c 

p  '•  dAo„  9  P 

A  °g  9 


(4.1-29) 


4.2  The  Intermediate  Functions 

The  mathematical  definitions  of  these  functions  are  as  follows  [subscripts  m 
and  m-1  denote  block  numbers]. 

The  grain  size  at  which  crack  initiation  occurs  is  calculated  from 


/ 


of'o' 


^  D  'm-r 


-1 


(4.2-1) 


which  is  written  by  solving  Eq.  3.2-la  for  D.  The  quantity  p(D)/0  is  independent  of  grain 
size  and  calculated  in  4.4.  Because  of  hardening,  D  is  not  always  the  largest  grain  size. 


aI(D,  Acp)  =  aI(  E  ,  ACp)  vs  0. 

^  m  PVl  ^ 


(4.2-2) 


I(D.AeJ)  »  1(0.  AeJ)  +  aI(D.  AeJ) 

^  m  ^  m-l  ^ 


m 


(4.2-3) 


I 


Aep"'*^(D)  is  the  solution  to 
m 
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1(0,  =  1  . 

m 

aH(D)„  *  aH(E„„  )  vs  0. 

Aa(D)^  *  1*  °* 

a(D)m  =  a(D)^_j  +  Aa(D) 


(4.2-4) 

(4.2-5) 

(4.2-6) 

(4.2-7) 


at  m  =  0,  a  =  a^.  A  good  value  for  a^  would  be  20  ym  for  the  A1  2219-T851  alloy  used  in 
our  study.  Initiation  requires  a(D)  «  a^,  where  a^  is  the  critical  initiation  crack  length 
which  is  equal  to  the  grain  depth. 

4.3  The  Cumulative  Functions 
These  are: 

H(D,ni)  =  (aH{D)  +  H(D,  tn-1))  vs  m  (4.3-1) 

Aep"'**(0,  m)  =  AeV’^(O)  VS  m  .  (4.3-2) 

4.4  Local  Mechanical  Properties 

All  bands  initiated  before  m*  blocks  where 


H(D,  m)  -  H(D,  m*)  =  1 


(4.4-1) 
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have  ceased  to  slip.  Wherefore 

=  Aep'""(D.  m*)  .  (4.4-2) 

AEp”*^'^(0)  and  AEp'"**(D)  allow  Epg(D)  p(0)  to  be  found.  Epg(D)  is  given  by  Eq. 
(3.2-6)  and 

Ae 

=  J  ^  n(A£p)  dAsp .  (4.4-3) 

*  ^min 

Acp 

4.5  Initiation  Lifetime 

In  its  present  form,  the  model  just  described  is  automatically  calibrated  for 
mean  and  sequence  effects  (for  A1  2219-T851)  by  the  microscopic  deformation  assess¬ 
ments.  While  all  of  the  lifetime  behavior  could  be  handled  this  way,  it  is  probably  best  to 
couple  some  of  the  model  calibration  to  actual  fatigue  data.  The  simplest  way  to  do  this 
is  to  exercise  the  model  for  constant  and  fully  reversed  stress  amplitude  and  for  several 
values  of  y  (in  Eq.  (4.1-24))  in  the  range  of  1-2  (reasonable  theoretical  estimates). 

A  comparison  of  predicted  and  experimentally  determined  initiation  lifetimes 
would  then  be  used  to  select  the  best  values  for  A  and  y.  Thereafter,  for  any  spectra 
applied  to  that  alloy,  A  and  y  used  in  lifetime  calculation  would  remain  the  same. 

The  model  will  predict  a(D).  This  actually  presupposes  that  there  will  be  at 
least  one  surface  grain  of  size  0  present  in  the  stressed  area,  containing  a  constituent 
particle  near  the  maximum  size.  Hence,  a(D)  «  gives  the  cycles  to  initiation.  If  a 
large  surface  area  is  stressed,  this  will  always  be  the  case.  But,  as  schematically 
illustrated  in  Fig.  4.5-1,  if  the  stressed  area  is  small  the  forecast  from  this  model  will  be 
conservative,  simply  because  the  sites  most  vulnerable  to  initiation  will  not  often  be 
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present.  In  this  case,  a  statistical  calculation  of  the  most  probable  initiation  lifetime 
could  be  made  and  would  require  knowledge  of  the  distribution  in  surface  grain  sizes  (Fig. 
4.5-2)  and  the  stressed  surface  area.  For  completeness,  the  model  must  include 
stochastic  variations  in  a^,  a^,  and  a^,  as  mentioned  earlier. 


•ean>« 


Fig.  4.5-1  Effect  of  surface  area  on  initiation  lifetime. 

4.6  Step-by-Step  Instructions  for  Modeling  Initiation  Lifetime  in  A1  2219-T851 

4.6.1  The  Cycle-by-Cycle  Calculations 

The  purpose  of  the  eye le-by -cycle  analysis  is  to  determine  the  local  damage 
rate  (over  one  loading  block)  due  both  to  stress  and  plastic  strain  ranges  for  a  grain  of 
arbitrary  size  and  mechanical  properties.  This  calculation  needs  to  be  done  for  just  one 
block,  and  allows  us  to  determine  the  incremental  change  in  mechanical  properties,  in 
hidden  damage  condition  (degree  of  accumulative  hardening  of  active  bands,  and 
predamage  leading  to  new  bands),  and  in  crack  extension. 

1.  Choose  from  the  flight  spectra  a  repetitive  load  block.  (This  will  look 
like  Fig.  4. 1 -3a,  but  for  the  applied  stress  rather  than  local  stress.) 
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Decode  this  into  a  sequence  of  minimum  and  maximum  stresses 

a  a  a  .a  a  _a 

^  mr  1^*  ^  fi 


_a  ^a  ^a  ^a  „a  ^a 

®j^min*  jmax*  ^2"*^* . * 


GRAIN  SIZE 


-EARLIEST  POSSIBLE 
INITIATION 


h— 10M— 
(D) 


•AFTER  10M  BLOCKS 


Fig.  4.5-2 


(a)  Evolution  of  the  a(0)  functions  with  cycles;  and  (b)  resulting 
probability  of  initiation  vs  grain  size. 
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For  each  half-cycle,  determine  the  local  stress/strains. 


3a.  Begin  in  the  grain  body  with  o  =  0.75,  o^J  =  140  MPa, 

n  9 


Eq'=  100,000  MPa,  Og  =  0,  ej  =  0  . 


3b.  Calculate  •  ®g  *  ®g  ^°g  half-cycle  using 

Eqs.  (4.1-1)  to  (4.1-10).  This  also  defines  a 


min  max 

sequence  ,  Og  . . . 


min  max 

M  f  * 

^n  ®n 


3c.  Consider  the  surface  layer  with  =  0.96,  E  =  40,000  MPa, 


Oj  =  70  MPa. 


pg 


3d.  Use  the  sequence  from  3b  to  find  Aep  half-cycle-by- 

half-cycle  using  Eqs.  (4.1-11)  to  (4.1-19).  For  each  half-cycle,  also 

find  Ojj,  Oj2»  calculation  on  the 

next  half-cycle.  This  will  give  sets  of  results  which  are  displayed 
similar  to  the  results  in  Fig.  4.1 -3a.  Both  stress  and  strain  range 
results  must  be  compiled  and  this  must  be  done  for  a  selected  set 
of  possible  values  of  Epg  in  the  grain's  main  body. 

3e.  Apply  "real-time"  rain  flow  algorithm  to  reduce  AOg,  Acp  ,  to 

equivalent  fully  reversed  spectra.  Truncate  and  Og  to  their 
positive  values  and  apply  rain  flow  to  convert  these  positive  going 
equivalents  as  well. 


3f.  At  the  completion  of  the  block,  arrange  the  elements  of  each  of 
the  five  spectra  in  ascending  amplitude  and  define  the  P 
distribution  function  and  the  intra-block  cycle  number  for  each,  as 
the  example  in  Fig.  4.1-3b. 
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3g.  Repeat  3b-3f  for  new  values  of  E  .  It  is  likely  that  results  for 
values  of  10'',  10^,  lO**,  10^  10  and  1  MPa  will  suffice. 

4.  Use  the  P's  to  find  AH,  Al  and  Aa  in  terms  of  Ep  using  Eqs.  (4.1-21), 
(4.1-23)  and  (4.1-29). 


4.6.2  The  Blodc-by-Blodc  Calculations 

The  block-by-block  calculations  deal  with  the  progressive  evolution  of  the  the 
damage  state  space,  and  are  supported  by  the  property  evolution  results  of  the  cycle-by¬ 
cycle  analysis.  The  spectra  response  is  embedded  in  the  cycle-by-cycle  analysis.  This  is 
just  like  a  simpler  conventional  summation  of  cumulative  damage  (i.e.,  1/n-type  damage) 
except  that  the  grain  size,  band  densities,  band  heU'dening  state,  and  predamage  to  form 
band  states  make  the  calculation  multidimensional. 

5.  Begin  with  the  first  block.  Set  a  =  a^,  E^'  =  100,000  MPa. 

5a.  Use  Eq.  (4.2-2)  to  redefine  aI  in  terms  of  grain  size  D.  [On  the 
first  pass,  4l(Epg,  Ae^)  is  independent  of  grain  size.] 

5b.  Find  I(D,  Acp)  and  solve  Eq.  (4.2-4)  for  Acp*"*^  .  For  each  grain 
size,  this  generates  results  as  in  Fig.  4.6-1,  defining  what  fraction 
of  potentially  viable  sites  for  b«md  formation  have  produced  bands. 

5c.  Likewise,  find  aH(D),  Aa(D)  and  a(D)  using  Eqs.  (4.2-4)  through 
(4.2-7). 

5d.  Compare  a(D)  to  a^.  to  test  for  initiation. 

5e.  Find  the  cumulative  functions  H  and  AEp"’®^(0)  [Eqs.  (4.3-1)  and 
(4.3-2)]j  and  use  these  and  Eqs.  (4.4-1)  and  (4.4-2)  to 
Ep  ^(0)  •  See  Fig.  4.6-2  for  the  construction  used  to  find 
for  one  grain  size. 
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5f.  Calculate  Ep^  and  p(D)/D  [Eqs.  (3.2-6)  and  (4.4-3)]  and  return  to  5a 
to  do  the  next  block. 

4.7  Model  Calibration 

The  bulk  of  the  material  coefficients  needed  to  calibrate  this  model  have  been 
defined  by  microscopic  deformation  observations  for  A1  2219-T851.  Now  that  these 
elements  have  been  identified,  only  a  few  weeks  laboratory  work  would  be  required  to 
obtain  the  calibration  values  for  a  new  alloy.  This  would  be  accomplished  by  producing  a 
sequence  of  measurable  strains  by  fatigue  (e.g.,  Ref.  5)  from  which  the  local  stresses, 
state  of  constraint,  and  hence  local  flow  stresses  and  strain  hardening  coefficient  can  be 
deduced.  Growth  parameters  y,  A  and  a^  must  be  selected  by  comparison  of  the  model 
predictions  to  measured  initiation  lifetimes  for  fully  reversed  loading  before  the 
spectrum  loading  calculations  can  be  done.  Likewise,  should  be  determined  by  a 
comparison  of  Ep^  to  experimentally  determined  values  vs  cycles  at  fixed  amplitude 
using  data  such  as  in  Section  2.1. 
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•C4«2n 


Schematic  illvjstration  of  procedure  used  to  find  the  range 

to  for  which  slip  bands  have  both 

initiated  and  have  not  yet  hardened.  This  defines  (with  n) 
the  total  number  of  active  slip  bands  at  block  m  in  a  grain 
of  size  D. 
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5.0  DISCUSSION 


Mean  stresses  might  affect  a  model  material  such  as  we  have  just  described  in 
several  ways  by: 


1 .  Altering  the  local  plastic  strain  ranges. 

2.  Changing  the  relationships  of  plastic  strain  to  the  response  functions  aH, 
Al  and  Aa.  The  first  avenue  is  directly  embedded  in  our  half-cycle  by 
half-cycle  calculation  of  the  local  deformations.  We  have  also  proposed 
a  truncation  of  Ae^  prior  to  calculating  Aa  in  an  attempt  to  care  for  a 
mode  I  nature  of  the  microcrack  growth.  Our  experiments  are  insuffi¬ 
cient  to  tell  if  AH  and  Aa  respond  to  mean  stress  in  addition  to  the  local 
plasticity.  However,  investigators  working  along  the  same  line  with  bulk 
low  cycle  fatigue  behavior  have  concluded  that  this  kind  of  mean  effect 
is  small,  if  it  is  present  at  all. 

Effects  of  load  sequence  in  our  model  will  appear  through 

1.  Their  impact  on  the  local  plastic  strains. 

2.  Because  of  an  interaction  with  the  hardening  state  H. 

The  extension  of  initiation  lifetime  of  samples  coaxed  by  low  amplitude 
loading  is  directly  attributable  to  the  accelerated  hardening  of  large  compared  to  small 

Q 

grains  in  the  alloy  surface.  While  a  conservative  lifetime  prediction  (and  simpler 
numerical  analysis)  might  result  if  hardening  is  omitted  from  the  model,  it  is  our  surmise 
that  an  important  potential  impact  of  load  sequence  would  also  be  lost  from  the 
predictions. 

To  what  kinds  of  materials  and  loading  situations  would  our  model  best  apply? 
First,  the  deformation  should  be  quasi-elastic  in  the  bulk  with  strains  less  than  10~^,  for 
most  of  the  load  sequence.  We  hope,  but  it  remains  to  be  demonstrated,  that  occasional 
plastic  overloads  will  not  compromise  our  predictions,  as  long  as  the  percentage  of 
lifetime  lost  in  the  extremes  is  small.  Second,  the  plasticity  should  be  localized  and. 
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hence,  constrained.  This  might,  for  instance,  happen  in  large  grains  in  high  strength 
structural  alloys.  It  might  also  happen  in  soft  spots  that  occur  in  some  superalloys  and 
titanium  alloys.  Structural  aluminum  alloys  in  a  T6  or  TS  heat  treatment  condition,  for 
instance,  should  be  representable.  A  large  grained  superalloy  like  WASPALOY  and 
certain  a-B  titanium  alloys  are  also  candidates.  Also,  a  fine-grained  material  such  as 
IN  100,  "if  it  has  soft  spots",  is  a  potential  candidate,  although  the  slip  band  portion  of  the 
model  will  probably  need  to  be  reworked  for  this  material.  If  the  alloy  is  cyclically 
stable  or  cyclically  hardens,  is  employed  below  its  bulk  yield  strength,  and  if  it  displays 
fatigue-induced  changes  in  surface  residual  stress,  indicative  of  localized  flow,  it  is  a 
candidate  for  this  analysis. 

What  about  possible  alterations  to  the  slip  band  nucleation  model?  Is  it  likely 
that  the  shallow  surface  layer  of  nucleation  sites  inferred  for  A1  2219-T851  is  ubiqu- 
tious?  Probably  not.  Our  surmise  is  the  confinement  of  band  nucleation  to  a  surface 
layer  occurs  only  in  coarse  precipitation  hardened  alloys  such  as  the  overaged  2219.  Yet 
to  be  established  by  experiment  is  the  prospect  that,  in  alloys  where  the  hardening  phase 
is  more  easily  cut  by  dislocations,  the  nucleation  sites  for  slip  bands  lie  more  deeply 
below  the  surface.  In  the  limit  that  band  nucleation  can  occur  within  the  entire  grain 
volume,  our  model  requires  a  few  mild  alterations,  which  will  actually  simplify  the 
numerical  calculations. 

1.  The  "surface  layer"  will  become  the  grain  volume  with  an  appropriate  a 
(<s  0.75  for  pancake-shaped  grains,  and  «  0.4  for  equiaxed  grains). 

2.  4l  will  be  calculated  for  the  grain  volume  using  the  "kinematic  model". 

3.  For  equiaxed  grains  Ep^  «  D~^  ,  as  the  entire  grain  contributes  to  band 
initiation. 

4.  The  local  stress  state  must  be  calculated  using  a  composite  model  of 
deformation  as  discussed  in  Section  2.2. 

How  might  the  crack  initiation  and  lifetime  controlling  factors  vary  with  alloy 
processing  and  with  location  in  a  structure.  We  believe  that  this  approach  will  handle 
those  alloys  where  soft  spots  in  the  surface  matter.  Materials  such  as  cast  iron,  which 
contain  initial  flaws  of  grain  size,  will  fail  by  crack  growth  at  rates  relatively  insensitive 
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to  local  plasticity.  But,  the  initiation  lifetime  of  many  of  the  high  strength  alloys  clearly 
depend  upon  the  micromechanics  of  local  deformation.  The  effect  of  grain  shape  (and 
surface  layers)  on  a  (the  proportiorwility  constant  between  plastic  strain  and  reaction 
stress)  and  hence  on  the  relative  importance  of  the  local  reaction  stresses  to  plastic  flow 
should  largely  determine  how  important  load  sequence  and  mean  stresses  will  be  to  life¬ 
time.  If  the  reaction  stresses  are  always  small,  because  grains  are  deep  and  wide  normal 
to  the  principal  stress  axis,  the  local  and  external  loads  will  be  nearly  equal  and  sequence 
should  be  less  important.  Of  course,  the  amplitude  content  of  the  spectra  and  its 
consequence  to  AH  and  Al  will  always  matter.  Reaction  stresses  will  normally  increase 
lifetime  by  lowering  the  plastic  strain  range.  So,  for  identical  spectra,  local  loading  of  a 
rolled  plate  in  the  transverse  direction  with  an  exposed  transverse  surface  will  likely 
accelerate  initiation  because  of  the  very  low  a  values  of  the  pancake-shaped  grains  for 
that  loading  orientation. 

Finally,  what  might  alloy  microstructure  do  to  the  crack  initiation  criterion 
and,  hence,  to  the  model  we  have  described?  Constituent  particles  common  to  structural 
aluminums  will,  if  sufficiently  large,  nucleate  cracks  which  progress  to  the  grain 
boundary.  If  the  alloy  is  overaged,  the  growth  is  more  likely  to  be  in  Mode  I  and  given  by 
our  model  exactly  as  it  is  described.  In  this  case,  the  largest  typical  particle  depth  would 
be  appropriate  to  use  as  the  starting  crack  size.  If  the  surface  area  under  stress  is  very 
small,  a  treatment  of  the  statistical  effects  of  particle  size  might  be  appropriate.  If  the 
aluminum  alloy  is  in  a  peak  hardness  heat  treatment,  crystallographic  growth  from  the 
initiation  site  is  more  likely.  should  not  be  truncated  before  Aa  is  calculated  if  this 
is  the  case  because  the  shear  bands  can  slip  in  compression  as  well  as  in  tension.  Growth 
to  the  boundary  of  a  crystallographic  crack  is  likely  to  be  mixed  mode  and  respond  to  the 
local  compressive  as  well  as  to  the  local  tensile  stress.  In  "clean"  alloys  having  only  very 
small  particles,  nucleation  and  growth  of  cracks  in  the  slip  bands  themselves  is  likely. 
Then,  we  suspect  that  a  microcrack  growth  analysis  based  on  the  quantity  Ae^/p  might 
be  most  appropriate,  as  this  is  proportional  to  the  strain  in  each  individual  band. 

Still  unresolved  is  how  to  handle  those  alloys  in  which  grain  boundary  cracks 
are  the  fatal  precursors.  In  A1  2219-T851,  isolated  boundary  cracks  form  before  the 
surface  softens,  apparently  in  response  to  the  microscopic  plasticity  of  the  surface 
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layer.  These  cracks  are  very  shallow  (a  few  microns)  and  have  not  been  found  to 
propagate  to  failure.  But,  the  mechanism  which  controls  this  behavior  remains  to  be 
identified. 
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6.0  SUMMARY 


In  many  structural  alloys,  crack  initiation  takes  place  at  the  surface  in  sites 
locally  softened  by  fatigue.  These  areas  behave  just  as  material  at  the  root  of  a  notch 
when  exposed  to  stresses  exceeding  yield.  The  plastic  flow  is  constrained  by  the  elastic 
member  surrounding  the  notch,  limiting  the  deformation  and  creating  reaction  stresses. 
These  can  cause  stresses  at  the  notch  root  to  differ  substantially  from  the  externally 
applied  stresses.  Load  sequence  and  mean  stress  then  alter  lifetime  because  the  reaction 
stress  modifies  the  coupling  between  load  sequence  and  local  plastic  strain  range. 

For  individual  softened  grains  in  an  alloy,  this  picture  is  further  complicated 
by  the  progressive  change  in  strain  softening/hardening  behavior  of  each  grain  at  rates 
dependent  on  grain  size.  In  A1  2219-T851,  one  result  is  that  cracking  is  focused  in  large 
grains  at  high  cyclic  stresses,  and  in  somewhat  smaller  (although  still  larger  than  the 
mean  size)  grains  at  lifetimes  close  to  the  fatigue  limit. 

While  we  have  a  reasoneible  insight  into  the  deformation  constraints,  there  are 
clearly  some  loose  ends  in  our  overall  understanding.  In  particular,  the  critical 
experiments  to  unequivocally  demonstrate  the  presence  of  a  thin  surface  layer  of  slip 
band  nucleation  in  A1  2219-T851  remain  to  be  done.  Of  more  importance  is  the  fact  that 
the  mechanics  of  the  property  evolution  in  single  grains  is  largely  understood  only 
phenomenologically.  The  science  which  might  allow  one  to  make  modifications  in  the 
model  of  property  evolution  to  recognize  failing  cases,  etc.,  remains  to  be  done. 

Nature  is  complex;  but  aircraft  design  must  go  on.  Consequently,  we  have 
taken  what  we  have  learned  about  local  surface  deformation  and  the  resulting  cracking 
and  embedded  this  in  an  initiation  lifetime  model.  Our  hope  is  that  when  this  is 
calibrated  using  a  minimum  of  test  data,  the  resulting  predictive  power  will  be  superior 
to  more  conventional  empirical  damage  models.  One  reason  for  this  hope  is  that  our 
focus  has  been  on  understanding  the  underlying  deformation  which  most  assuredly 
controls  initiation.  Our  understanding  of  the  interplay  between  creation  and  hardening  of 
slip  bands  as  molded  by  the  deformation  constraints  is  also  fairly  comprehensive. 
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Just  as  with  contemporary  models,  we  propose  the  use  of  rain  flow  counting  to 
reduce  the  calculated  local  stress  and  plastic  strain  spectra  to  a  collection  of  fully 
reversed  cycles  for  a  loading  block.  But,  this  must  be  done  to  three  different  spectra, 
namely,  the  local  stress  and  strain  in  the  grain  body  and  plastic  strain  in  the  surface 
layer  within  each  grain.  Each  spectra  is  determined  by  a  simple  half-cycle-by-half -cycle 
calculation  of  stresses  and  strains  within  a  grain  over  the  duration  of  the  block.  The 
representation  of  these  results  in  terms  of  a  range  of  possible  plastic  moduli  during 
individual  blocks  give  a  set  of  response  functions  describing  increments  in  damage 
involving  slip  band  nucleation,  slip  band  hardening  and  incremental  microcrack  growth 
for  a  loading  block.  These  functions,  then,  can  be  applied  in  sequence  to  predict  the 
initiation  lifetime  to  formation  of  a  grain-sized  crack  vs  the  local  grain  size. 

Mathematical  models  of  the  local  grain  stiffness,  and  the  constrained  plastic 
strain  range  are  stated,  and  a  phenomenological  description  of  the  slip  band  nucleation 
and  band  hardening  rates  vs  the  local  plastic  strain  range  or  stress  range  are  stated. 
Data  to  allow  these  submodels  to  be  calibrated  for  A1  2219-T851  are  presented.  The 
piece  that  is  missing  from  this  work  is  the  numerical  framework  to  generate  the  loading 
block  response  functions.  The  equations  to  assemble  these  functions  are  stated  and 
algorithms  to  do  a  "real-time"  rain  flow  processing  of  the  local  spectra  are  available. 
Only  the  numerical  means  to  construct  smooth  functions  which  include  local  stiffness 
variables  in  the  three  response  functions  must  be  formulated  to  create  a  working  model. 
Thus,  all  the  essential  model  elements  have  been  defined,  but  the  computer  code  to 
implement  the  initiation  lifetime  calculation  for  an  arbitrary  spectra  must  still  be 
written. 
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7.0  RECOMMENDATIONS 

This  model  is  ready  to  be  encoded  and  tested  with  a  comparison  of  predictioru 
to  actual  initiation  lifetime  vs  loading  spectra.  If  the  model  performs  reasonably,  it 
should  be  formulated  into  a  user-friendly  design  tool.  The  next  steps  are: 


1.  To  find  out  how  well  the  model  works  for  several  classes  of  structural 
materials.  The  important  variables  and  considerations  are: 

Grain  size 

Crain  aspect  ratio 

Surface  layer  presence/absence 

Deformation  mode  (slip  bands  vs  local  Mode  I  cracking) 

Stress  amplitude  content 
Load  sequence 
Mean  stress 

Statistical  (and  notch)  effects  with  an  appropriate  model  refinement 

2.  To  identify  appropriate  ways  to  calibrate  the  model  by  determination  of 
responses  to  standard  load  spectra,  in  order  to  minimize  the  information 
needed  in  the  local  deformation  behavior  and  grain  size  distribution. 

3.  To  explore  alternate  mathematics  to  implement  this  model.  We  have 
described  a  quasi-sequential  calculation  of  average  lifetime.  Other 
avenues  possibly  computationally  more  efficient  are  state  space 
evolution  and  probabilistic  formulations. 
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APPENDIX  A 


Appendix  A  was  presented  at  the  Fatigue  S7  Conference  in  Charlottesville, 
VA,  June,  19S7.^  The  notation  is  essentially  the  same  as  through  the  main  body  of  this 
report.  The  figures  and  references  for  this  section  are  self-contained. 
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THE  EVOLUTION  OF  LOCAL  MECHANICAL  PROPERTIES  OF  A1  2219-T851 
DURING  FATIGUE 


W. L.  Morris  ,  M.R  James*  and  B.N.  Cox 


Local  strain  amplitudes  were  measured  In  Individual 
grains  of  an  A1  2219-T851  alloy  during  fatigue. 

The  alloy  has  a  360  MPa  bulk  cyclic  yield  strength, 
but  Its  surface  Is  microplastic  at  cyclic  stresses 
greater  than  70  MPa.  With  fatigue  at  ±  275  MPa,  an 
upper  flow  stress  of  -  200  MPa  develops  In  grains 
much  larger  than  the  mean  size.  Large  residual 
stresses  caused  by  constraint  of  the  localized 
deformation  severely  limit  the  total  strain  at 
these  softened  sites.  The  use  of  a  load  reduction 
sequence  to  minimize  these  residual  stresses  during 
strain  measurements,  so  as  to  obtain  more  accurate 
values  for  the  local  flow  stresses.  Is  discussed. 


INTRODUCTION 


Large  local  residual  stresses  are  always  present  during  Inho¬ 
mogeneous  deformation  (1,2).  In  fatigue-softened  alloys,  loca¬ 
lized  deformation  Is  modified  by  residual  stresses  resulting  from 
constraint  of  flow  by  the  harder  surroundings  (3,4).  Even  on  a 
smooth  surface,  deformation  Inhomogeneity  at  the  grain  level 
ensures  the  presence  of  these  residuals,  and  hence  a  substantial 
difference  between  the  local  and  externally  applied  stresses.  A 
simple  analysis  has  shown  (3)  that  the  residual  stress  In  Indi¬ 
vidual  grains  can  easily  be  so  large  that  portions  of  a  grain  can 
be  In  compression  while  the  external  stress  Is  tensile.  This 
complicates  a  simple  determination  of  the  local  flow  stress  of 
Individual  grains  In  alloys. 

In  two  recent  */aper»  (3,4),  we  used  models  of  the  deformation 
of  a  soft  ellipsoid  In  an  elastic  matrix  to  calculate  the  flow 
stress  and  strain  hardening  characteristics  of  surface  grains  In 
A1  2219-T851.  The  models  account  for  the  local  residual  stress, 
allowing  for  the  Interpretation  of  strains  measured  In  individual 


Published  in  Fatigue  87,  ed.  R.O.  Ritchie  and  E.A.  Starke, 
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grains  as  a  function  of  external  stress.  We  find  that  fatigue  of 
A1  2219-T851  at  stresses  below  Us  bulk  cyclic  yield  strength 
causes  a  reduction  In  the  local  0.2X  offset  yield  strength  of  the 
interiors  of  the  grains,  with  preference  for  the  largest  surface 
grains.  Grains  3-8  times  larger  than  the  60  wb  aean  are  the  most 
affected  and  are  also  the  eventual  sites  of  crack  Initiation. 
Ultimately,  the  yield  strength  of  these  grains  falls  to  half  the 
bulk  cyclic  value,  and  well  below  the  applied  cyclic  stress 
amplitude.  However,  the  local  strain  amplitude  typically  remains 
less  than  10'^  in  these  softened  grains  because  the  deformation  Is 
constrained  by  the  nearly  elastic  surface. 

These  results  seem  In  conflict  with  ample  evidence  that  many 
alloys  locally  harden  (5,6).  Even  In  A1  2219-T851,  Knoop  indenta¬ 
tion  shows  that  fatigue  hardens  the  Interiors  of  Individual 
grains,  with  the  change  being  greatest  In  the  largest  grains  (Fig, 
1).  Previous  local  strain  measurements  In  this  alloy  have  shown 
that.  In  large  grains,  the  width  of  the  local  strain-external 
stress  hysteresis  loop  at  zero  stress  reaches  a  peak  with  fatigue 
and  then  decreases  in  amplitude,  again  suggestive  of  hardening 
(7). 

In  this  paper,  we  demonstrate  that  the  small  local  flow 
stresses  deduced  from  strains  measured  within  Individual  grains 
are  not  an  artifact  of  our  deformation  analysis.  While  strain 
hardening  parameters  and  deformation  depth  cannot  be  accurately 
obtained  without  a  rigorous  theoretical  analysis,  the  local  flow 
behavior  can  be  found  directly.  A  load  reduction  sequence  is  used 
to  minimize  the  local  residual  stress,  causing  the  external  and 
local  stresses  to  be  nearly  equal.  This  technique  Is  used  to 
characterize  the  evolution  with  fatigue  of  flow  stresses  and 
strain  hardening  at  low  plastic  strain  amplitudes  In  a  300  urn 
grain  In  A1  2219-T851.  The  apparent  conflict  between  Indentation 
and  local  strain  measurements  Is  not  completely  resolved,  but  the 
new  local  stress-strain  data  provide  some  clues  to  possible 
answers  to  this  problem.  We  believe  that  the  10%  strains  encoun¬ 
tered  during  Indentation  (8)  can  make  this  method  a  misleading 
probe  of  the  highly  constrained  yielding  which  occurs  In  isolated 
softened  grains.  In  which  the  maximum  plastic  strains  are  only 
10”5  during  typical  fatigue  experiments. 


EXPERIMENTAL  TECHNIQUES 


The  360  MPa  yield  strength  A1  2219-T851  had  a  60  vm  average  grain 
size  In  the  rolling  plane  measured  at  45*  to  the  rolling  direc¬ 
tion.  The  grains  vrere  pancake  shaped,  being  nominally  20  urn 
deep.  Tapered  flexural  specimens  (9)  of  the  alloy  were  carefully 
machined  and  polished  to  minimize  Initial  residual  surface 
stresses  and  then  chemically  etched  to  reveal  the  grain  boun- 
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darles.  Fatigue  was  In  bending  at  nominally  t  275  MPa  In  dry 
nitrogen.  The  local  strains  reported  were  measured  using  a 
reference  gauge  placed  within  the  center  of  grains  and  at  least 
25  ym  from  the  nearest  boundary  by  a  micromanipulator,  and  with 
the  aid  of  a  loading  jig  In  a  scanning  electron  microscope  (SEM), 
as  described  In  Ref.  10.  The  reference  gauge  Is  a  thin  mica  flake 
used  as  a  convenient  ruler  and  attached  to  the  substrate  electro¬ 
statically  so  that  the  surface  can  deform  freely  beneath  It.  It 
acts  as  a  stable  length,  needed  to  compensate  for  the  ±  5  >  10'> 
magnification  Instability  of  our  SEN.  Distortion  and  parallax 
errors  from  the  ruler  are  smallest  at  zero  load,  where  we  made 
most  of  our  measurements.  For  fully  reversed  loading,  the  width 
of  the  local  strain-external  stress  hysteresis  loop  at  zero  load 
was  found  by  comparing  a  micrograph  taken  after  completing  the 
compressive  cycle  to  one  taken  after  the  next  tensile  cycle.  The 
Information  needed  to  deduce  flow  stress,  strain  hardening  and 
constitutive  relationships  was  obtained  by  changing  the  loading 
sequence  prior  to  the  loop  width  measurement  (4).  We  used  a 
stereoscopic  analysis  of  the  micrographs  to  maximize  the  sensitiv¬ 
ity  of  measurements  of  the  surface  displacements  relative  to  the 
ruler.  Since  the  displacements  of  nearly  Identical  high  contrast 
objects  in  each  pair  of  micrographs  were  found,  the  accuracy  was 
much  better  than  the  point-to-point  resolution  of  the  SEM.  The 
displacement  of  any  high  contrast  profile  can  easily  be  determined 
to  a  small  fraction  of  the  profile  width.  It  Is  the  Image  repro¬ 
ducibility  which  actually  matters.  The  quality  of  our  SEM  Images 
has  been  progressively  improving  over  the  past  few  years  and  we 
can  now  achieve  about  ±  15A  in-plane  displacement  sensitivity. 


RESULTS 


To  determine  how  the  flow  stress  of  a  grain's  Interior  progres¬ 
sively  changes  during  fatigue,  we  Interrupted  the  constant 
amplitude  cycling  at  Intervals,  and  determined  the  local  flow 
stress  by  measuring  strains  within  the  grain  for  a  series  of 
reduced  cyclic  loads.  The  local  mechanical  properties  In 
A1  2219-T851  change  slowly,  so  this  measurement  activity  has  no 
significant  effect  on  the  property  evolution.  The  cyclic  stress 
was  dropped  to  a  value  o  and  six  cycles  were  applied  to  allow  the 
local  stress-strain  hysteresis  to  reach  equilibrium  (4).  Then  the 
residual  strain  at  zero  external  load  was  measured  over  a  tensile 
loading  cycle.  The  strain  found  at  each  o  is  the  width,  W^,  of  an 
external  stress-local  strain  hysteresis  loop  at  zero  exter^l 
load.  will  always  be  smaller  than  the  zero  load  width  of  a 
local  stress-local  strain  loop  because  of  local  residual  stress. 

As  the  measurement  stress  o  Is  decreased,  U.  at  equilibrium 
becomes  progressively  smaller.  At  the  same  time,  the  maximum 
local  residual  stress  within  the  grain  decreases  because  the 
maximum  local  plastic  strain  amplitude  has  decreased.  For  a 
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sufficiently  small  plastic  strain  amplitude,  the  difference 
between  the  external  and  local  stresses  at  equilibrium  will  be 
small,  and  an  accurate  determination  of  the  local  flow  stress  can 
be  made  from  W.  vs  o.  An  upper  bound  on  the  local  residual  stress 
magnitude  is  given  by  the  product  of  the  local  plastic  strain  and 
the  alloy  modulus  (3).  For  aluminum,  the  modulus  is  nearly  inde¬ 
pendent  of  grain  orientation,  and  the  maximum  plastic  strain  is 
rarely  5  times  larger  than  W^;  so  for  W^^'s  smaller  than  5  »  10"*, 
the  difference  between  external  and  local  stress  should  be  less 
than  20  MPa  after  equilibrium  Is  reached.  Only  resolved  shear 
stress  variability  due  to  grain  orientation  should  remain. 

High  sensitivity  measurements  of  made  early  in  fatigue  show 
the  surface  is  microplastic  at  stresses  greater  than  70  MPa  (Fig, 
2).  In  fact,  the  stress  dependence  of  for  a  very  large 
(300  urn)  and  a  smaller  (120  urn)  grain  are  essentially  identical. 

We  would  expect  to  see  more  scatter  in  such  results  from  resolved 
shear  stress  with  more  measurements  but,  clearly,  the  deformation 
in  the  essentially  unfatigued  alloy  is  short  ranged  in  that  it 
does  not  sense  the  grain  boundaries. 

With  continued  fatigue,  a  dramatic  change  in  loop  width  occurs 
in  just  the  largest  grains.  Figure  3  compares  the  behavior  of  a 
300  viti  grain  after  100  and  250  cycles.  Within  the  t  1.5  «  10* ^ 
measurement  sensitivity,  there  is  little  change  in  the  lower  flow 
stress  (oj,  at  Wq  •  1  »  10*i)  with  fatigue.  But,  characteristic 
of  large  grains  in  the  alloy  is  the  appearance  of  an  upper  flow 
stress,  Oy,  in  this  case  seen  clearly  after  250  cycles.  Below 
the  grain  has  apparently  strain  hardened  with  fatigue  (i.e., 
dw^/do  has  decreased);  above  Oy,  the  grain  has  strain  softened. 

We  have  shown  the  results  in  Fig.  3  on  a  logarithmic  scale  to 
clarify  the  data  at  small  strains.  When  these  and  similar  data  at 
other  fatigue  increments  are  plotted  linearly,  trends  in  o,  and  c^j 
with  fatigue  are  more  accurately  defined  (Fig.  4a),  and  indirect 
information  on  the  strain  hardening  behavior  can  also  be  obtained 
(Fig.  4b)  by  treating  the  strain  as  linear  with  o  above  and  below 
Oy.  A  theoretical  analysis  of  loop  widths  for  such  two-stage 
yielding  is  presented  in  Ref.  4. 

We  see  from  Fig.  4a  that  the  local  flow  stresses  are  nearly 
constant  in  a  300  vm  grain  with  perhaps  a  small  fatigue  induced 
decrease  in  Oy.  The  important  changes  with  fatigue  are  in 
dWy/do.  The  quantity  dW^/do  Is  shown  for  the  stationary  yield 
surface  approximation  In  Ref.  4  to  be  equal  to  2(l-o)/(2Ep  *  aE.); 
where  a  has  a  value  less  than  1  and  Is  a  proportionality  between 
the  local  residual  stress  and  the  plastic  strain  determined  by 
deformation  depth,  grain  shape,  and  Poisson's  ratio.  E.  and  E. 
are  plastic  and  elastic  moduli.  An  Increase  in  dW^/do  means  that 
either  Ep  or  a  has  decreased.  Typical  values  for  a  will  be 
between  i  and  0.6,  with  very  shallow  deformation  raising  o.  A 
full  interpretation  of  the  dW^/do  results  in  Fig.  4b  requires 


A5 

C8577D/bje 


NADC-89044-60 

FATIGUE  87 


Rockwell  International 

Science  Center 


SC5470.fr 


additional 


Wq  vs 


load  sequence  data  so  that  a  and  E.  can  be 


Independently  determined  (see  Ref.  4).  However,  If  the  depth  of 
deformation  within  the  grain  is  reasonably  constant,  the  results 
mean  that  between  100  and  500  cycles  Ep  has  Increased  below  and 
has  decreased  above  Og,  i.e.,  the  grain  has  strain  hardened  below 
o„  and  strain  softened  above  a^..  A  more  complete  analysis,  to  be 
described  later,  suggests  that  the  fall  In  dWp/do  at  800  cycles 
and  above  Og  Is  a  consequence  of  a  reduction  in  the  depth  of  the 
deformation  within  the  grain.  The  results  shown  at  500  cycles  are 
for  a  different  300  ym  grain,  so  considerably  more  data  are  needed 
to  completely  sort  out  the  fat1gue>1nduced  evolution  of  local 
strain  hardening  and  deformation  depth. 


DISCUSSION 


It  is  simplistic  to  think  of  the  fatigue-induced  changes  in  local 
surface  deformation  of  A1  2219  as  reflecting  traditional  hardening 
or  softening.  The  changes  are  complex  and  partially  hidden  by  the 
large  residual  stresses  developed  by  constraint  of  the  local 
deformation.  Analysis  in  Ref.  4  shows  that  by  500  cycles  Og  is 
just  below  the  conventional  (0.2*  offset  strain)  yield  strength. 
The  loop  width  measured  at  maximum  cyclic  stress  is  small  because 
the  deformation  is  constrained.  Mughrabi's  (11)  concept  of 
treating  the  flow  of  such  a  material  as  that  of  a  composite  of 
matrix  and  penetrating  channels  or  slipbands  seems  relevant 
here.  But  unlike  the  case  of  Cu,  it  appears  to  us  that  the  upper 
rather  than  the  lower  flow  stress  more  likely  corresponds  to  long- 
range  dislocation  motion.  The  argument  for  this  Is  Indirect. 
First,  Og  develops  after  fatigue.  While  PSB's  are  not  optically 
visible  in  A1  2219-T851,  banded  slip  has  been  found  by  high 
resolution  strain  field  analysis  in  large  grains  after  fatigue 
(3),  suggestive  of  a  connection  to  Og.  Second,  before  there  is 
substantial  fatigue  damage,  the  flow  above  Oj  is  independent  of 
grain  size  and,  therefore,  probably  short  range  and  in  the 
matrix. 


No  substantial  change  in  the  flow  characteristics  of  small 
grains  has  been  found  after  10-  cycles  at  ♦  275  MPa,  while  In 
300  ym  grains  changes  both  above  and  below  Og  are  apparent  after  a 
few  hundred  cycles.  So,  If  there  Is  strain  hardening  below  Og  In 
the  matrix,  it  must  have  accompanied  the  formation  of  the  long- 
range  dislocation  structure  which  deforms  above  Og,  and  cannot 
Instead  be  the  cause  of  the  long-range  structure.  Otherwise, 
matrix  strain  hardening  would  take  place  over  the  entire  surface 
early  In  the  lifetime. 


How  might  an  Indentation  test  perceive  a  composite  grain  In 
which  the  matrix  had  strain  hardened,  and  In  which  banded  slip  had 
developed  for  which  flow  above  Og  up  to  strains  of  perhaps  10"^ 
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was  nearly  perfectly  plastic?  Clearly,  on  the  basis  of  a  lower 
0.2%  offset  yield  strength,  fatigue  has  softened  the  larger 
grains,  but  there  are  several  ways  that  Indentation  night  sense 
hardening.  The  strains  at  the  tip  of  a  pyramidal  Indenter  are 
relatively  Insensitive  to  indentation  depth  and  are  nearly  10~‘ 
(8).  Any  mechanism  which  exposes  the  matrix  to  Indentation  and 
suppresses  the  contribution  of  flow  from  the  bands  will  sense 
fatigue  hardening.  Possibilities  for  this  Include  a  saturation  of 
the  strain  In  the  bands  at  the  much  higher  strain  of  Indentation, 
or  the  destruction  of  the  long-range  slip  mode  by  the  Indenta¬ 
tion.  Naturally,  we  do  not  believe  that  softening  is  being 
universally  misinterpreted  as  hardening.  We  have  recently  found 
unequivocal  local  cyclic  hardening  of  an  A1-4%  Cu  alloy  with  the 
same  local  strain  measurement  techniques.  Our  results  simply  call 
for  care  in  the  analysis  of  local  mechanical  properties. 

We  have  known  since  the  work  reported  in  Refs.  3  end  4  that 
deformation  in  large  (300  pm)  grains  in  A1  2219-T851  eventually 
becomes  nearly  perfectly  plastic  during  fatigue  (i.e.,  E-  «  0) . 
While  it  has  been  clear  that  reaching  this  state  involves  a  dra¬ 
matic  drop  in  the  local  yield  strength,  it  was  not  clear  if  the 
evolution  of  the  local  mechanical  properties  entailed  a  decrease 
in  Oy,  a  decrease  in  or  changes  in  both  these  quantities.  It 
now  appears  likely  that,  while  both  may  be  changing,  the  most 
important  effect  is  on  E  . 


SUMMARY 


Strains  measured  in  individual  grains  of  A1  2219-T851  show  a 
progressive  evolution  of  the  local  mechanical  properties.  The 
entire  surface  is  microplastic  at  fully  reversed  stresses  greater 
than  ±  70  MPa,  Early  in  fatigue,  the  deformation  strains  are 
independent  of  grain  size,  suggesting  that  the  flow  resides  in  a 
matrix  subject  to  short-range  constraints.  With  continued 
fatigue,  an  upper  flow  stress  Oy  becomes  visible  at  -  200  MPa  in 
grains  substantially  larger  than  the  mean  size.  Below  Oy,  the 
material  apparently  strain  hardens,  above  Oy  it  strain  softens. 
The  possibility  that  the  development  of  Oy  coincides  with  the 
penetration  of  the  matrix  by  slip  bands  is  suggested. 
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CRAIN  SIZE,  A<m 

Figure  1  Knoop  hardness  of  A1  2219-T851  vs  the  Individual  grain 
size  In  a  single  specimen  shows  apparent  hardening  after  fatigue. 
Each  datum  is  an  average  of  approximately  ten  measurements. 


CYCLIC  (TMtS  Mr* 

Figure  2  Comparison  of  loop  widths  measured  at  zero  load  vs 
cyclic  stress  amplitude  for  a  120  and  a  300  ym  grain.  The  data 
were  obtained  before  a  substantial  fatigue  Induced  change  In 
deformation  has  occurred. 
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Figure  3  The  progressive  change  in  loop  width  in  a  300  um  grain 
with  fatigue. 
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Figure  4  Trends  in  (a)  Oj  and  o^;  (b)  dW^j/do  above  and  below 
for  the  300  um  grain. 
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APPENDIX  B 


Appendix  B  was  presented  at  the  Spring  1987  meeting  of  the  TMS-AIME  in 
Denver,  CO,  in  February,  1987.^  The  notation  used  is  essentially  the  same  as  in  the  main 
body  of  this  report.  The  figures  and  references  in  this  section  are  self-contained. 
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LOAD  SEQUENCE  EFFECTS  ON  THE  DEFORMATION  OF  ISOLATED 


MICROPLASTIC  GRAINS 


M.R.  James  and  W.L.  Morris 

Rockwell  International  Science  Center 
1049  Camino  Dos  Rios 
Thousand  Oaks,  CA  91360 


Abstract 

Strains  measured  in  individual  large  grains  of  A1  2219-TS51  are  used  to  deduce  the 
local  constitutive  behavior  that  controls  plastic  flow  under  spectrum  loads.  A 
composite  model  for  the  deforming  grain  is  most  consistent  with  experiment.  The 
material  apparently  contains  two  intermingled  dislocation  structures  whose  flow  is 
controlled  by  different  yield  criteria.  Special  load  sequences  are  used  to  determine  the 
yield  criterion  best  descriptive  of  each  component,  based  on  strains  measured  over 
100  um  gauge  lengths  in  individual  grains. 
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Even  for  "elastic"  loading,  the  role  of  inhomogeneous  localized  plasticity  in 
fatigue  has  been  demonstrated  repeatedly  (1-3).  Recent  theories  of  inhomogeneous 
plasticity  in  tensile  and  cyclic  loading  principally  consider  large  strains  and  show  the 
importance  to  plastic  flow  of  inhomogeneity  in  alloy  microstructure  (4,5)  (e.g.,  granu¬ 
larity)  and  of  local  instabilities  in  deformation  (6).  When  dislocation  phenomena  have 
been  used  to  model  the  constitutive  response  of  local  elements  in  a  material,  it  is  often 
supposed  that  the  underlying  events  of  dislocation  generation,  saturation  and  recovery 
are  sequential  and  conditiorted  on  the  local  total  strain  range.  However,  for  the  small 
plastic  strains  common  to  most  structural  fatigue  failires,  the  dislocation  structures 
themselves  tend  to  be  inhomogeneous.  Dislocation  channels  or  bands  penetrate  a 
matrix  of  dislocation  tangles,  creating  a  microscopic  composite  (7-9).  The  local  stress 
state  in  such  subelements,  their  local  mechanical  properties  and  constitutive  response 
may  all  differ,  and  each  quasi-homogeneous  material  element  will  "remember”  prior 
deformation  and  react  to  factors  beyond  total  local  strain  and  stress.  A  complex 
interaction  between  load  s^uence,  local  stress-strain  response  and  rates  of  fatigue 
crack  initiation  can  be  anticipated  for  all  such  materials. 

Recent  studies  of  the  deformation  of  individual  large  grains  in  an  overaged  alumi¬ 
num  alloy  (10-14)  suggest  the  importance  of  a  "composite"  dislocation  structure  in 
modeling  the  local  strain  -  external  stress  response.  The  flow  of  each  isolated 
microplastic  grain  is  constrained  by  the  nearly  elastic  surface,  creating  reaction 
stresses  reminiscent  of  those  at  the  root  of  a  locally  plastic  notch.  Procedures  have 
been  developed  (10,11)  to  determine  the  mechanical  properties  of  isolated  plastic  grains 
from  strains  measured  over  microscopic  «  100  vm)  gauge  lengths.  The  local  internal 
stresses  that  might  affect  the  measured  property  values  are  handled  by  load  shedding, 
which  either  reduces  the  stresses  (12),  or  allows  the  resulting  plastic  hysteresis  to  be 
analyzed  (10,11).  After  fatigue  of  A)  2219-T851,  marked  hardening  of  the  plastic  flo* 
is  observed  below  a  critical  stress,  above  which  there  is  simultaneous  and  dramatic 
strain  softening  (12).  Does  this  softening  herald  the  development  of  banded  slip,  seen 
only  by  high  resolution  strain  field  mapping  (13)  in  2219,  or  is  it  better  explained  by  the 
saturation  and  recovery  mechanism  of  a  single  subcomponent  of  the  material,  utilized 
by  Estrin  and  Kubin  (6)  in  modeling  inhomogeneous  plasticity? 

To  resolve  this  question  we  examine  the  local  constitutive  behavior  of  individual 
grains  by  studying  their  response  to  several  load  sequences.  In  particular,  we  determine 
what  constitutive  behavior,  in  terms  o(  motion  of  the  yield  locus,  must  be  assigned  to 
each  element  of  a  two-component  composite  in  order  to  achieve  the  observed  stress- 
strain  reaction  to  mean  stress  and  tensile  reloadings  found  experimentally. 


Backaround 

When  Al  2219-T851  is  loaded  cyclically  below  its  360  MPa  cyclic  yield  strength, 
the  initial  strains  in  a  300  vm  grain  found  over  a  tensile  increment  are  small,  but 
discernible  above  the  experimental  measurement  error,  at  stresses  greater  than  70  MPa 
(o.)  (Fig.  la).  A  few  hundred  cycles  at  t  270  MPa  changes  the  local  properties  of  the 
grain,  creating  en  upper  flow  stress  (o„)  near  200  MPa,  strain  hardening  below  o^,  and 
softening  above  (rig.  2a)  (12).  The  te^nique  used  to  measure  strains  over  the  luO  um 
gauge  lengths  required  for  this  study  (see  Refs.  10,  11,  15,  16)  employs  displacements 
measured  between  two  points  of  zero  external  load  in  a  loading  sequence.  The  total 
strains  calculated  from  the  displacements  reflect  reactions  to  local  plasticity  (residual), 
rather  than  the  elastic  strains  stemming  from  an  external  load.  If  the  residual  strains 
are  less  than  -  4  >  10**  the  reaction  stresses  will  typically  be  sufficiently  small  that 
data  such  as  in  Fig.  2a  will  depict  the  true  local  stress-strain  response.  If  the  residual 
strains  are  large,  deformation  models  are  needed  to  separate  the  elastic  reaction  strain 
and  plastic  strain  components,  in  order  to  identify  flow  points  and  calculate  strain 
hardening  coefficients. 
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Other  studies  on  A1  2219-T83J  (10,11)  indicate  that  in  the  overaged  state,  local 
mechanical  properties  do  not  vary  greatly  across  a  grain  so  that  an  assumption  of 
uniform  plastic  deformation  is  adequate  for  fully  reversed  loading.  Two  limiting  cases 
of  the  flow  surface  locus  have  been  examined  for  microplastic  grains,  namely:  1)  a 
"stationary"  surface  (S)  independent  of  prior  load;  and  2)  a  ^nematic"  surface  (K), 
which  moves  freely  with  the  maximum  stress  excursion  on  the  previous  load  reversal 
(11).  With  the  S  model,  the  grain  has  no  memory  of  beir^  deformed  other  than  its 
instantaneous  plastic  strain  and  the  resulting  reaction  to  its  constraint  by  the  grain 
boundaries,  while  with  the  K  model  a  maximum  possible  ^dc  stress"  is  remembered  in 
addition  to  plastic  strain  and  boundary  constraints.  The  deformation  of  microplastic 
grains  for  which  and  a|  are  both  sutionary  or  both  kinematic  has  been  considered 
elsewhere  (11).  The  presence  of  a  K  modelled  behavior  for  has  been  ruled  out  as 
giving  an  unacceptable  description  of  strain  transients  following  a  sudden  change  in 
cyclic  stress  (11),  but  the  possibility  of  a  K  model  of  and  an  S  model  of  has  not 
been  examined. 


For  this  study  we  model  the  latter  situation  by  defining  a  composite  material  for 
the  grain,  assumed  to  experience  the  same  local  stress  in  each  element.  Each  of  the 
two  components  (i  s  1,2)  is  assigned  a  "plastic”  modulus  for  plastic  strain  above  its 
respective  flow  stress,  giving  a  plastic  strain  increment 

(1) 

The  strains  in  the  two  elements  are  summed;  thus  E-  contains  a  hidden  dependence  on 
both  the  strain  hardening  and  the  relative  volume  of  ^e  two  components.  The  flow  cri¬ 
teria  are  applied  independently  to  each  element  in  a  numerical  calculation  (based  on 
equations  in  Ref.  10)  of  the  strains  and  stresses  at  each  load  reversal.  This  allows  o.  to 
be  defined  according  to  a  K  model  and  to  an  S  model  status.  (The  complete 
formulation  of  th.s  model  v,ili  be  discussed  elsewhere.) 

In  addition  we  have  considered  a  mode)  describing  a  nonuniformly  deforming  grain 
represented  by  an  interior  enveloped  by  a  shell  of  different  mechanical  properties  near 
the  gram  boundary.  The  stress-strain  behavior  of  such  dual-domain  systems  has  been 
described  by  Morris  et  al.  (10).  The  constitutive  equations  for  a  dual-domain  grain  viere 
implemented  numerically  for  an  arbitrary  load  sequence.  In  the  following  we  use  U  and 
0  to  denote  uniform  and  dual-domain  gram  models;  the  entry  K/S  refers  to  a  K  model 
of  Oj  and  an  S  model  of  o  .  For  example,  D-S/S  is  a  dual-domain  model  with  a 
stationary  flow  surface  locus  Tor  both  and  o^. 


Results 


Deformation  in  a  300  um  grain  before  develops  during  fatigue  is  examined  in 
Fig.  1.  The  sample  was  cycled  at  t270  MPa  for  100  cycles  and  then  the  load  was 
dropped  to  e.  as  shown  in  the  insert  of  Fig.  la.  Measurement  of  was  made  after  3 
cycles  at  to  allow  for  equilibrium  of  the  local  stress-strain  response.  Repetition  of 
load  shedding  for  multiple  values  of  has  been  found  most  useful  to  determine  the 
local  mechanical  properties  of  the  grain  (11).  Material  values  were  chosen  for  U-S/S 
and  U-K/S  models  to  fit  the  observed  rcsiAial  strains  in  Fig.  la.  However,  for  these 
same  parameters,  only  the  U-K/S  model  depicts  the  strain  behavior  (or  a  sequence  of 
tensile  reloads  to  successively  higher  tensile  maxima  (o^iax^ 
increase  in  the  (low  locus  that  occurs  with  a  K  model  is  necessary  to  prevent  signifi¬ 
cant  grain  deformation  on  reloading. 

After  a  few  hundred  fatigue  cycles,  Oy  appears  (Fig.  2a).  Residual  strains 
measured  during  load  shedding  for  fully  reversed  loading  are  again  described  by 
selection  of  appropriate  material  parameters  for  the  two  regions  of  flow  (or  the  U-K/S 
and  U-S/S  models.  Using  these  coefficients,  only  the  U-K/S  model  correctly  represents 
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the  increa'.'ng  strain  in  response  to  a  series  of  stress  relaxation  cycles  around  a' 
increasingh  more  positive  mean  stress  <Omean^  ^  kinematic  increase  of  the 

0|  flow  locus  suppresses  the  large  internal  strains  that  would  develop  with  increasing 
mean  stress.  Another  way  this  result  might  occur  is  if  flow  at  the  grain  boundary 
lowered  the  maximum  stress  experienced  by  the  grain  interior.  This  possibility  was 
examined  using  the  dual-domain  model  (10).  With  appropriate  selection  of  material 


properties  (or  the  boundary,  the  D-S/S  model  can  almost  mimic  the  behavior  in  Fig.  2b 
but,  for  the  same  material  parameters,  the  representation  of  the  behavior  in  fully 
reversed  loading  (Fig.  2a)  is  poor. 


Fig.  1  Local  strain  behavior  in  a  3CD  urn  gram  after  100  cycles  at  t  270  MPa.  VI  is  a 
residual  strain  measured  at  zero  externa!  load  in  the  load  sequence  uolid 
circle),  relative  to  an  earlier  point  (open  circle),  a)  As  a  function  of  the 
applied  cyclic  stress  amplitude;  b)  for  a  series  of  loads  of  increasing  tensile 
maximum  following  unloading  from  270  MPa. 


Fig.  2  Local  strain  behavior  in  a  300  um  grain  after  230  cycles  at  *  270  MPa.  is  a 
residual  strain  measured  at  zero  external  load  (solid  circle)  relative  to  an 
earlier  point  (open  circle),  a)  As  a  function  of  the  applied  cyclic  stress 
amplitude;  b)  following  a  stress  envelope  converging  to  a  positive  mean  stress 
"mean*  to  an  envelope  collapsed  to  a  zero  mean. 
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Discussion 

In  a  grain  of  monolithic  properties,  if  is  less  than  for  any  loading  sequence, 
then  the  relationship  must  hold  for  all  sequences.  This  is  simply  because  each  region  of 
flow  in  a  monolithic  material  will  be  conditioned  by  particular  levels  of  strain.  To  en¬ 
sure  this  behavior  in  a  two-stage  yielding  model,  both  the  lower  and  upper  flow  stresses 
would  have  to  be  either  stationary  or  kinematic.  In  contrast,  the  strain  behavior  of 
single  grains  we  report  for  the  mean  stress  and  tensile  reloading  experiments  is  repre¬ 
sented  well  only  if,  (1)  the  high  strain  hardenii^  deformation  found  just  above  occurs 
at  a  stress  less  than  in  fully  reversed  loading,  and  (2)  if  is  greater  than  e^  after 
single  large  tensile  excursions.  It  appears  that  for  this  to  happen,  must  be  kine¬ 
matic,  o^  stationary,  and,  further,  that  the  grain  must  consist  of  a  microscopic  compos¬ 
ite  of  material  elements  of  two  types,  one  having  a  flow  stress  the  other  o^.  It  also 
suffices  to  treat  the  deformation  as  uniform  throughout  the  grain,  so  the  individual 
composite  subcomponents  must  be  small  compared  to  a  100  um  gauge  length.  Since  a 
stationary  flow  locus  is  tantamount  to  a  large  separation  between  dislocation  pinning 
points,  is  Oy  then  associated  with  banded  slip?  This  would  mean  that  onset  of  flow  in 
the  bands  ^n  require  higher  stresses  than  flow  in  the  matrix  (above  o^),  an  uncom¬ 
fortable  proposition  except  that  flow  in  the  bands  involves  much  larger  strains  than  in 
the  matrix. 

For  Ai  2219-T851  it  is  becoming  clear  that  the  relationship  of  to  e^  in  indi¬ 
vidual  grains  is  extremely  dependent  on  the  prior  load.  It  is,  therefore,  unavoidable 
that  the  local  plastic  strains  in  this  material  will  be  highly  sensitive  to  the  details  of 
load  sequence  as  well  as  to  the  amplitude  range.  U  ith  the  assistance  of  a  simple  model 
depicting  the  way  the  material  properties  in  the  grain  evolve  with  cumulative  plasticity 
(1<*),  it  should  now  be  possible  to  calculate  how  load  sequence  will  effect  the  local 
strains,  and  hence  the  progress  of  fatigue  crack  initiation  in  isolated  microplastic 
grains. 


Summary 


The  complexities  of  the  deformation  of  a  large  microplastic  grain  in  Al  2219-T851 
are  revealed  by  local  strain  measurements  made  for  external  loading  other  than  fully 
reversed.  Two  stages  of  flow  are  recognized,  and  a  comparison  with  models  of  the  de¬ 
formation  process  suggests  the  likelihood  that  each  takes  place  in  a  distinct  subcompo¬ 
nent  of  the  grain.  Many  of  these  elements  must  be  intermingled  to  form  a  microscopic 
composite,  as  the  deformation  appears  to  be  quasi-uniform  over  gauge  lengths  compar¬ 
able  to  the  grain  size.  Oie  subcomponent,  known  to  strain  harden  quickly  with  fatigue, 
is  best  described  by  a  kinematic  movement  of  the  yield  surface  locus  with  prior  loading, 
suggesting  a  reaction  to  short  range  back  stresses  in  the  grain.  Deformation  of  the 
material  is  therefore  extremely  dependent  on  the  prior  load  sequence.  The  second  com¬ 
ponent  dramatically  strain  softens  with  fatigue,  and  is  best  described  by  a  "stationary” 
yield  surface  locus  -  suggestive  of  a  long-ranged  constraint  of  flow  by  grain  bound¬ 
aries.  A  physical  representation  of  our  model  equations  might  be  a  matrix  containing 
slip  bands,  although  we  know  for  the  2219  aluminum  alloy  that  traditional  persistent  slip 
bands  are  not  present.  Nevertheless,  the  single  uniform  deformation  model  which 
incorporates  both  of  these  kinds  of  behavior  provides  an  adequate  description  of  the 
deformation  of  the  isolated  nonplastic  grains  in  the  2219  aluminum  alloy  that  are 
responsible  for  fatigue  crack  initiation. 
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